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ABSTRACT 


Two  methods  of  accounting  for  body-lif ting  surface  interference  in 
unsteady  flow  are  considered.  The  first  method  is  described  in  Part  I 
of  this  report,  while  the  second  will  be  described  in  Part  II  to  follow. 

The  first  method  is  a  direct  application  of  nonplanar  lifting  surface 
elements  to  both  the  lifting  surfaces  and  the  body  surfaces.  The  body  is 
treated  as  an  annular  wing.  This  type  of  idealization  must  be  used  with 
an  axial  doublet  introduced  to  account  for  body  incidence  effects.  The 
undesirable  effects  of  the  annular  wing  representation  are  then  reduced. 

The  second  approach,  to  be  described  in  Part  II,  uses  an  image  system 
and  an  axial  singularity  system  to  account  for  the  effects  of  the  bodies. 

This  report  also  describes  an  improvement  of  the  Doublet-Lattice  Method 
of  Albano  and  Rodden.  The  improvement  pertains  to  wing-tail  problems  where 
there  exists  a  small  vertical  (non-zero)  separation  between  the  wing  and 
tail  planes.  Such  problems  can  now  be  handled  with  ease. 

Tni.j  volume  contains  the  development  of  the  theory,  correlation  of 
theory  with  experimental  data,  and  the  parametric  study.  Volume  II  is  a 
guide  to  the  computer  program  and  contains  the  FORTRAN  listing. 
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1.0  INTRODUCTION 


This  report  presents  a  direct  application  of  the  Doublet-Lattice  Method  to 

problems  involving  body-lifting  surface  interference.  The  Doublet-Lattice  Method 

of  Albano  and  Rodden*  which  is  used  here, is  an  extension  of  the  steady  Vortex - 

2  3  4 

Lattice  Method.  Developments  by  Hedman,  Giesing,  and  James  have  shown 
that  the  Vortex-Lattice  Method  is  (1)  simple  (2)  accurate  and  (3)  quite  versatile. 
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Correlations  by  Kalman,  Rodden  and  Giesing  show  that  the  same  attributes 
apply  to  the  unsteady  Doublet-Lattice  Method.  A  survey  of  the  development  of 
the  steady  and  unsteady  lattice  methods  is  given  in  Reference  5. 

Two  methods  of  accounting  for  body-lifting  surface  interference  in  unsteady 
flow  are  considered.  The  first  method  is  described  in  Part  I  of  this  report, 
while  the  second  will  be  described  in  Part  II  to  follow. 

The  first  method  is  a  direct  application  of  nonplanar  lifting  surface  elements 
to  both  the  lifting  surfaces  and  the  body  surfaces.  The  body  is  treated  as  an 
annular  wing.  Such  an  idealization  has  been  used  by  Woodward^  in  the  steady 
case.  This  type  of  idealization  must  be  used  with  an  axial  doublet  introduced  to 
account  for  body  incidence  effects.  The  undesirable  effects  of  the  annular 
wing  representation  are  then  reduced. 

The  second  approach,  to  be  described  in  Part  II,  uses  an  image  system  and 
an  axial  singularity  system  to  account  for  the  effects  of  the  bodies.  Basically  it 
will  be  an  extension  of  the  method  of  Reference  3. 

This  report  also  describes  an  improvement  of  the  Doublet-Lattice  Method 
of  Albano  and  Rodden.  The  improvement  pertains  to  wing-tail  problems  where 
there  exists  a  small  vertical  (non-zero)  separation  between  the  wing  and  tail 
planes.  Such  problems  can  now  be  handled  with  ease. 


2.  0  THEORETICAL  DEVELOPMENT 


2.  1  The  Doublet-Lattice  Method 


The  velocity  normal  to  an  oscillating  surface,  W  =  UmRe 


related  to  the  lifting  pressure  Ap  =  q  Re(ACp  e*wt),  by  the  integral  equation 


/  iwt\  . 
\w  e  ),  i  s 


w(x,  y,  z) 


=  Jj'  K(x-|,  y-n,  z-;,o,  M)  ACp  d£  dr 
L.S, 


(2.1-1) 


where  £  is  the  streamwise  coordinate,  cr  is  the  tangential  spanwise 
coordinate  (see  Sketch  2.  1-1),  w  is  the  frequency  of  oscillation,  M  is  the 
Mach  Number  and  K  is  the  kernel.  The  limit  L.  S.  indicates  integration  over 
all  the  lifting  surfaces.  The  integral  in  Eq.  (2.  1-1)  can  be  approximated 
as  follows: 


w(x,  y,  z) 


_1_ 
8  ir 


E  ACp 

s  b 


K(x-£,  y-n,  z-4,w,  M)  d£  d  cr 


(2.  1-2) 


Element 

s 

where  s  is  an  index  indicating  the  sending  element  A  further  approximation  is 
introduced  in  the  Vortex-  and  Doublet-Lattice  Methods.  The  integration  of  K 
in  the  streamwise  direction  (  £)  is  done  simply  by  lumping  the  effect  into 


SKETCH  2.  1-1  WING  AND  ELEMENT  COORDINATES 


a  loaded  line  at  the  1 /4-chord  line  of  the  element. 


w(x,  y,  z)  =  £  AC  —  —  f  K(x- y-1,  z-C,  w,  M)dcr  (2.1- 

5  Ps  Sir  J  1/4 

Element 

s 


3) 


The  result  is  an  unsteady  horseshoe  vortex  whose  bound  portion  lies  along 
the  1  /4-chord  line  of  the  element. 


In  this  equation  the  normalwash  boundary  condition  w(x,  y,  z)  is  known 

and  the  lifting  pressure,  AC  ,  over  each  element  is  unknown.  A  set  of 

P  s 

linear  algebraic  equations  may  be  formed  from  (2.  1-3)  if  the  normalwash  is 
satisified  at  as  many  points  as  there  are  elements.  Equation  (2.  1-3)  can 
then  be  written  in  matrix  form 


w 


[D]  *Cj 


(2.  1-4) 


where  a  typical  element  of  [D],  D^,  is 


D 


rs 


/ 


K(x-Si/4,  y-n,  z  -i,  w,M)do- 


(2.  1-5) 


Element 


s 


and  where  x,  y,  z  (both  for  Dfs  and  w^)  are  the  coordinates  of  the  receiving 
point  "r".  There  is  one  control  or  receiving  point  per  element  and  the  sur¬ 
face  normalwash  boundary  condition  is  satisfied  at  each  of  these  points.  The 
control  point  is  centered  spanwise  on  the  three-quarter  chord  line  of  the 
element  (Sketch  2.  1-2).  This  choice  of  control  point  location,  shown  by 

4 

James  to  be  optimum  for  the  two-dimensional  case,  results  in  a  high  degree 
of  accuracy.  Sketch  2.  1-2  also  shows  the  idealization  of  the  surface  into 
elements  or  boxes  which  are  arranged  in  strips  parallel  to  the  free  stream 
so  that  surface  edges,  fold  lines,  and  hinge  lines  lie  on  element  or  box  boundaries. 


OOWNW ASH 
COLLOCATION  POINT ■ 


Sketch  2.  1-2.  Idealization  of  a  Wing  Panel  into  Boxes 
2.  1.  1  Integration  Over  an  Element 


The  general  form  of  the  nonplanar  kernel,  K,  used  is  given  by  Rodemich 

g 

and  Landahl  as: 


K  -  p"1  r T" 

K.  —  e  Uco 


(KjT!  +  K2T2)/r' 


(2.  1-6) 


where  K, ,  T.,  K„,  and  T-,  are  given  in  Appendix  A.  The  basic  idea  of  the 
Doublet-Lattice  Method  is  to  fit  the  numerator  of  K  with  a  parabola  and 
integrate  (2.  1-5)  analytically.  Since  the  steady  kernel  can  be  integrated 
exactly  without  the  parabolic  fit  it  is  apparent  that  the  accuracy  will  be 
improved  if  only  the  unsteady  increment  is  integrated  using  the  approximating 


parabola.  Equation  (2,  1-6)  becomes 


K(x-|i/4>  y-n,  z-C,w,  M)  = 


,"1  Uco  (x'^l/4) 


e  Jc0  1/4  K1-K1(s) 


+  T, 


u 


e'iu”  (X"ei'4)K2-K2(s' 


+  TlKl(s)  +  T2K2 


(S)J/r2 


(2.  1-7) 


The  superscript  "s"  indicates  the  steady  case  (  w  =  0).  The  first  of  the 
bracketed  terms  in  the  numerator  of  Equation  (2,  1-7)  is  fitted  with  a  parabola 
and  integrated  while  the  second  of  the  bracketed  terms  is  integrated  analytically 
without  the  use  of  a  curve  fit. 


D 


rs 


D  (U  +  D  ^ 

rs  rs 


(2.  1-8) 


where 


D 


rs 


(1)  AX, 

8ir 


/' 


--2 

Ail  +  B  il  +  C 


dq 


D 


(s)  _  ^ 


rs 


(S)  +  T,K?(s) 

2  2  2 
I  r 


and  where  Aq2  +  Bq  +  C  =  Tx  [e1  “  (x"  h  /4}  -K^]  +  T2  [e1  uJX‘  h  /4}K2-K2(s)] 


Here  "e"  is  the  semi-width  of  the  element  in  the  plane  of  the  element.  The 
integration  involved  in  is  shown  in  Appendix  B.  The  term  D rJS^ 

represents  a  steady  horseshoe  vortex  and  its  evaluation  is  given  in  Appendix  C 
The  change  in  variable  from  <r  to  q  indicates  that  element  coordinates  are  to 
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be  used.  (See  Sketch  2.1-1)  Element  coordinates  for  a  particular  element 
are  centered  on  the  sending  element  and  rotated  into  the  plane  of  the  element. 
The  relations  between  element  coordinates  and  the  usual  coordinates  are 
given  as: 

x  =  x  -£c 

y  =  (y-nc)cosVs  +  (z-  t,c)  sin  Yg 

z  =  (z-  4C)  cos  Yg  -  (y-r)c)sinYs 

~r  =  Yr  "  Ys 

where  £  ,  U  and  £  are  the  coordinates  of  the  center  of  the  1 /4-chord  line 
c  c  c 

of  the  sending  element  and  Yg  is  the  dihedral  angle  of  the  element.  The 
expressions  for  and  £>  are  the  same  as  above  with  x,  y  and  z  replaced 

with  *1  and  4  . 

The  integration  used  to  obtain  D  ^  works  very  well  for  all  cases, 
planar  and  nonplanar,  with  one  exception.  Consider  the  case  where  the 
receiving  point  is  downstream  of  the  sending  element.  When  there  exists  a 
small  vertical  separation  between  the  receiving  point  and  sending  element 
then  the  numerator  of  the  kernel  has  large  variations  with  *1  and  a  parabolic 
fit  will  not  give  accurate  results. 

The  problem  can  best  be  understood  by  considering  the  nature  of 
the  kernel.  It  can  be  shown  that  the  kernel  is  a  semi-infinite  line  doublet 
whose  strength  varies  like  exp  [-iw(£-£j /^)/Uco3  t^‘e  point  This 

semi-infinite  line  doublet  is  to  be  integrated  in  the  ^-direction  (spanwise 
in  the  plane  of  the  element).  When  the  receiving  point  is  downstream  of  the 
sending  point  of  the  kernel  with  a  small  vertical  and  lateral  gap,  i.  e.  , 
small  "z  and  "y,  then  the  flow  field  is  dominated  by  the  local  strength  of  the 
semi-infinite  doublet  line. 
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Expressions  for  the  various  parts  of  the  kernel  valid  for  small  lateral 
distances,  and  downstream  of  the  sending  point,  in  element  coordinates, 
are: 


lim 

p2r/(*-61/4)  -  0 


x  >  6 1/4 

when 


K,  2 


K- 


-4 


=  cos  Y 


x2 


=  |z^/r^jcosYr  -|z(y-q)/r^jsin  Yf 


-  —  2  -? 

(y  -h)  +  z 


(2.  1-9) 


These  terms  may  be  used  to  construct  the  unsteady  increment  given  byr  the 
first  term  in  brackets  on  the  right-hand-side  of  Equation  (2,  1-7).  The 
numerator  of  this  term  is  to  be  fit  with  a  parabola.  Investigation  of  the 
variation  of  this  term  with  q  will  show  why  a  parabola  is  inadequate  and 
it  will  also  suggest  a  way  to  remedy  the  situation. 


Using  Equation  (2.  1-9)  above  in  the  numerator  gives: 


lim 

>2 


„(s)1.„.  r.'1  U»(x'^l/4)„  „  (s) 


K2-KHr 


_2 

z 


H2cos  Y  -  4  cos  Y  — T  + 
r  r  r2 


+  4  sin  "y  | 


(2.  1-10) 


The  terms  that  cause  the  trouble  arise  from  the  nonplanar  term  and  are 

2  _2  2  _ 2  — 

the  ones  that  are  divided  by  r  .  Plots  of  z  /r  and  -z(y-q)/r  (when  y 

is  zero)  are  given  in  Sketch  (2.  1-3). 
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SKETCH  2.  1  -3 


When  the  element  length,  2e,  Is  large  compared  with  z  the  variations  across 
the  element  are  very  large  and  a  parabola  is  inadequate. 

The  obvious  solution  to  this  problem  is  to  consider  the  nonplanar  terms 
separately  from  the  planar.  The  nonplanar  terms  have  basically  a  different 
behavior  than  the  planar  terms.  The  planar  terms  vary  like  1/r  but  the 
nonplanar  terms  vary  like  1/r4.  This  suggests  that  the  kernel  should  be 

v/ritten  as: 


K 


(Kj  Tj/r 


+  K2T|/r4) 


(2.  1-11) 


where 


-r* 

l2 


=  r  T 


=  [(z-Ucos  Y  '  (y-HJsin  Yr]  [(z-^)cos  Yg  -  (y-n)  sinYg) 
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Again  the  steady  kernel  should  be  considered  separately  to  increase  the 


accuracy* 


K  J(T1[e-1^6,'4,K1-K1«“']/r2  ♦  T»[e-iu»(X'6l^K2-K2<s»]/r4S 

+  jTjK^'/r2  +  T*  K2<S»/r4J  (2.  1-12) 


of  course,  the  steady  kernel  is  integrated  separately  to  a  steady  horseshoe 
vortex.  The  numerators  of  each  of  the  terms  in  the  first  bracket  may  now 
each  be  fit  with  a  parabola  since  they  are  very  slowly  varying  functions  of 


*1  across  an  element. 


n2  +  Bxn  + 


[•“ 


U«(3i'^l /4}k  .K2(s)]«A2n  +  Bzn  +  C2 


(2.  1-13) 


Equation  (2.  1-5)  then  becomes: 


D  =  D  (U  +  D  (2)  +  D  (s) 

rs  rs  rs  is 


(2. 1-14) 


where 


(2)  _ 


Al7'2  *  B2n  *  C1 


6  A2n2  +  b2t  +  c2 


2  T  3 
4 


=  f  (T1K1(8,/r2  +  T*K2(8,/r4)  dq 


The  integrals  D  ^  and  D  ^  are  the  same  as  those  encountered  in  Equation 
(2.  1-8)  except  that  the  coefficients  A,  B,  C  in  (2.  1-8)  are  not  the  same  as 
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(2) 

A,,  B,,and  C.  in  (2,  1-13).  The  term  D  '  '  is  new  and  represents  the 
ill  rs 

nonplanar  unsteady  increment.  The  evaluation  of  D  v  is  given  in 

r  s 

Appendix  II.  The  final  expression  is: 

(2)  (e/2)  A  |  (  2(y2+i2+e2)  (A2e2  +  C2)  + 4y  e2 

^rs  ~  -_2,_2  2.  )  r/-  .2  ,  -2i  r.-  .2  -2-, 

(y  +z  -e  )  (  [(y+e)  +  z  J  [(y-e)  +  z  ] 

-  "J  [(yVz2)A2+yB  2  +  C2]|  (2.1-15) 

where 

fL+ligj)  tan"1!  ifR . y  J 

2e|z|  (  (y  +z  -e  )  )  . 

Appendix  B  gives  an  expansion  for  a  valid  for  small  values  of  z  and  an 

-2  _2  2 

equation  valid  for  small  values  of  y  +  z  -  e  . 

2.1.2  Reflection  Planes  and  Ground  Effect 

Most  applications  of  the  present  method  involve  configurations  with 
one  or  two  planes  of  symmetry.  The  X-Z  plane  is  a  plane  of  symmetry 
for  most  cases,  e.  g.  ,  the  right  side  of  a  lifting  surface  is  the  mirror 
image  of  the  left  side.  In  other  cases  an  additional  plane  of  symmetry 
(the  X-Y  plane)  exists.  For  instance  a  lifting  surface  attached  to  a  ground 
plane  (as  in  a  wind  tunnel)  or  in  the  proximity  of  the  ground  (aircraft  in 
ground  effect)  represents  a  configuration  with  two  planes  of  symmetry. 

All  flow  conditions  can  be  split  up  into  symmetrical  and/or  anti- 
symmetrical  parts  in  relation  to  either  plane  of  symmetry.  If  the 
configuration  is  symmetrical  and  the  flow  is  either  symmetrical  or 
antisymmetrical  then  a  considerable  savings  in  computational  effort  can 
be  realized.  Consider  the  general  case  of  two  planes  of  symmetry.  The 
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entire  lifting  surface  system  is  made  up  of  segments  in  the  upper  right 
(UR)  quadrant,  the  upper  left  (UL),  the  lower  right  (LR)  and  the  lower 
left  (LL).  If  the  distribution  of  lifting  pressure  and  normalwash  are 
each  split  into  their  component  segments  (one  per  quadrant)  and  so  designated 
by  the  subscripts  UR,  TJL,  LR  and  LL  the  matrix  equation  given  by  Eq.  2.  1  -4 
can  be  written  as; 


dll  dl  2  d13  d14 

d2l  d22  d23  d24 

d31  d32  d33  d34 

d32  d33  d34  d35 

where  d..  are  submatrices  of  D... 

U  U 

Let  the  flow  symmetry  betweer  »*ight  and  left  be  indicated  by  the  term 

{1  symmetry 
-1  antisymmetry 
V  0  no  symmetry 

In  addition  let  «  indicate  symmetry  between  upper  and  lower. 

11  symmetry 

-1  antisymmetry  (ground  effect) 

0  no  symmetry 
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Then  AC  TTT 
pUL 

*  6iCpUR 

ACpLR 

=  £ACr>UR 

(2.  1-17) 

ACPLL 

*  6ACpLR 

= 

6<ACpUR 

The  distributions  for  all  of  the  quadrants  have  been  reduced  to  one 
distribution  for  the  upper  right  quadrant  Using  (2.  1-17)  in  (2,  1  - 1 6 ) 

gives 

WUR  =  [  dll  +  6d12  +  * dl 3  +  €6dl4  J  |ACpUr| 

Thus  only  that  part  of  the  configuration  that  appears  in  the  upper  right  hand 
quadrant  needs  to  be  considered. 

Configurations  such  as  T-tails  present  no  special  problems;  however, 
the  results  obtained  by  the  present  method  for  these  cases  must  be  properly 
interpreted.  The  vertical  surface  of  the  T-tail  is  in  reality  two  surfaces. 
Consider  the  configuration  in  Sketch  2.  1-4. 


Sketch  2.  1-4.  Symmetry  considerations  for  T-tails 
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The  T-tail  is  formed  in  the  limit  as  the  dihedral  of  surface  "a"  (Y  ) 

cl 

approaches  90°.  The  vertical  fin  is  then  made  up  of  surface  "a"  plus  its 
reflection  a'  .  The  present  method  determines  the  loading  only  on  sur¬ 
faces  in  the  right  hand  sector  and  thus  only  half  of  the  loading  is  determined 
for  the  vertical  fin.  The  lifting  pressure  distribution  and  the  spanwise 
loading  and  moment  are  off  by  a  factor  of  two.  (Program  H7WC  tests  for  this 
condition  and  adjusts  these  distributions.  See  Vol  II)  The  image  lifting  surface 
is  always  accounted  for  when  calculating  stability  derivatives  or  generalized 
forces.  Since  the  vertical  fin  surface  "a"  has  an  image  a1  the  derivatives  are 
calculated  correctly.  When  dealing  with  aeroelastic  problems  the 
structure  is  always  cut  in  half  along  a  symmetry  plane.  This  cut  splits  the 
vertical  fin  of  a  T-tail  into  two  equal  surfaces.  Only  half  of  the  aero¬ 
dynamic  forces  are  applied  to  each  half  of  the  fin.  Thus  the  AIC  calculation, 
to  be  described  in  Section  2.  5  for  use  in  aeroelctstic  problems,  is  formed 
properly  by  the  present  method. 

2.  2  Lifting-Surface/Body  Interference 

The  treatment  of  lifting-surface/body  combinations  presented  in  this 

£ 

report  (Part  I)  is  similar  to  that  of  Woodward  .  Woodward  considered 
only  the  steady  case  for  a  single  body  whereas  the  more  general  oscillatory 
case  for  multiple  bodies  is  dealt  with  here.  The  effects  of  the  isolated 
bodies  are  obtained  using  slender  body  theory.  Body -lifting  surface  inter¬ 
ference  is  obtained  by  placing  lifting  surface  elements  on  the  body  surface 
near  the  body-lifting  surface  intersection. 

The  order  of  solution  again  is  similar  to  that  employed  by  Woodward. 

The  slender  body  solution  for  the  longitudinal  force  distribution  on  all 

o 

bodies  is  obtained  first.  The  slender  body  solution  used  (J.  W.  Miles  ) 
handles  bodies  of  circular  cross-section  whose  radius  varies  in  the  longi¬ 
tudinal  direction.  The  motion  of  the  body  may  be  arbitrary.  The  longitudinal 
force  distribution  on  the  bodies  causes  a  flow  field  which  affects  the  lifting 
surface  elements.  The  normalwash  due  to  this  flow  field  is  then  determined 
at  all  lifting  surface  elements  except  those  on  the  body  surfaces.  The  axial 
singularities  found  for  a  single  isolated  body  using  the  slender  body  theory 
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divert  the  flow  around  that  body.  The  resulting  values  of  normalwash  found 
at  the  interference  lifting  surface  elements,  for  that  body,  are  therefore  zero. 
For  the  case  of  several  bodies  there  is  a  small  nonzero  normalwash  caused 
by  the  singularity  systems  of  the  other  bodies.  Exact  solutions  to  the  unsteady 
compressible  wave  equation  are  used  to  determine  the  flow  fields  generated 
by  the  bodies.  Specifically>pressure  doublet  distributions  are  used.  The 
normalwash  caused  by  the  body  axial  force  distributions  is  added  to  the  normal- 
wash  generated  by  the  motions  of  the  lifting  surface.  The  resulting  boundary 
value  problem  is  then  solved  exactly  as  outlined  in  the  last  Section  (2.  1). 

2.2,1  Slender  Body  Force  Distributions  and  their  Resulting  Flow  Fields 

q 

The  slender  body  theory  of  J.  W.  Miles  is  used  for  the  unsteady  case. 
The  longitudinal  lift  distribution  for  harmonic  motion  is  given  by  Equation 
16  of  Reference  9  and  is 


3L 

0  X 


m  (x )  w 


(f) 


iwt 

e 


(2.2-1) 


where  m(x)  is  the  virtual  mass  of  the  cross  section  and  w(x)  is  the  normal¬ 
ized  upwash"  (W^/Uco)  in  the  direction  of  L.  For  a  circular  cross-section 

2 

of  radius  Ro(x),  m(x)  =  l4pxRo  (x).  Introduction  of  this  value  of  m(x)  into 
Equation  (2.2-1)  gives: 


1  JLL  _L 

q  8x  2R 
^  o 


w<f>  +  R 


i 

w  + 


wR 


(2.  2-2) 


The  term  on  the  left-hand-side  of  Equation  (2.  2-2)  may  be  regarded  as  a 

lifting  pressure  coefficient,  AC  acting  on  the  projected  body  area.  The 

P 

force  per  unit  length  is  obtained  simply  as  this  pressure  multiplied  by  the 
local  body  diameter  and  dynamic  pressure. 


dL_ 

9x 


q  2R  AC  (f) 

n  o  p 


(2.  2-3) 


/ 
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and  thus  from  (2.2-2): 


AC 


(0  _ 


=  2i 


!  R  w(f>  +  w(fi'R  n  t  ik  w<f> 

|  o  or 


R  /c  1 
o 


(2.2-4) 


CO  C 


where  the  reduced  frequency  k  =  “rrr —  has  been  introduced. 

r  u  cn 


The  pressure  AC  actually  acts  only  on  the  body  axis  and  thus 
P  (f) 

acts  like  a  delta  function  whose  integrated  value  is  2  R  AC  .  Further- 

m  P 

more,  AC  '  has  two  components  in  general,  a  vertical,  aC7,  and  a 
P  (f)  L‘ 

horizontal,  ACy»  as  does  the  upwash  to  the  fuselage  w'  .  Thus 


Tc(f) 


6(y-\.  z-Sa) 


ACzk 


+  AC, 


where  rj  and  4  are  the  coordinates  of  the  body  axis.  Also  it  can  be  shown  that 
a  a 

if  a  pressure  doublet  is  placed  within  a  circular  body  the  total  integrated  force 
is  0.  5  (AC^^^AA);  i.  e.  ,  the  pressure  doublet  is  only  half  as  effective  when  it 
lies  within  a  circular  body.  Thus,  if  the  lifting  surface  theory  is  to  be  used  to 
obtain  the  flow  field  due  to  AC^  then  2AC^f^  must  be  used.  A  slight  generali¬ 
zation  of  (2.  1-1)  gives: 


wB(x,  y,  z)  = 


+ 


JL 

4t r 


1 

4tt 


// 

// 


Kz(x-£,  y-T  7.-1,,  w,  M) 


Ky(x-£,  y-T  z-4,u>,  M) 


6(y-na>  z_  £a)  ACZ  dtdr) 

(2.2-4) 

6(y-na,  z-^a)  ACy  d£d£ 


where  the  subscripts  Zand  Y  on  K  indicate  the  direction  of  the  pressure 
doublet  sheet,  and  the  B  subscript  on  w  indicates  that  the  normal  wash  is 
caused  by  the  bodies.  Integration  of  the  delta  function,  6,  in  the  1  and  £ 
directions  gives 


wB(x>  y.  z) 


+ 


(x-£,  y-  na,  z -  £a,  w,  M) 

Kv(x-£,  y-n  .  z-C  ,  w>  M) 
x  a  a 


AC,  2R  d£ 
Z  o 

AC Y  2Rq  d£ 


(2.  2-5) 
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Integration  in  the  longitudinal  direction  is  done  numerically.  Each  body 

axis  is  divided  into  a  number  of  elements  over  which  R  C,,  and  R  C 

o  Y  o  Z 

are  assumed  constant. 

wB(x,y,z)  =  2^  E  ^Aczs  Rog  J  Kz(x-|,  Y~\<  Z-Ca>«.  M)  d^j 


Element 


(2.  2-6) 


A  E  (aCy  R0g  /  Ky(x-^,  y-na,  M)d$) 


Element 


Integration  over  each  element  is  done  by  lumping  the  effect  at  the  leading 
edge  of  each  element,  i.  e.  ,  at  4  =  £ag 


j(x>y-z)  =  T^^'ACZ  Ro  H  KZ  (x"^a  *  y-na>  M) 


s  s  s  s  s 


(2.2-7) 


+  -zi  E  ACy  RQ  A|  Ky  (x-ea  ,y-ria,  Z-^.w.M) 


s  s  s  s  s 


The  values  of  the  normalwash  at  wing  (or  tail,  etc.  )  lifting  surface  control 
points  and  other  body  lifting  surface  control  points  may  be  assembled  into 
a  column  matrix  jwg  j  .  Using  (2.  2-7)  this  column  matrix  may  be  written  as 


{  wb}  =  tFz)  {Acz}  +  [fy1  {acy} 


(2.  2-8) 


where  typical  elements  of  the  non-square  matrices  [F„  ]  and  [Fv  ]  are 


tt-  R  Ae  K 
2tr  o  S-  Z 

s  s  rs 


(2.  2-9) 


27  R0  A^s  Ky 


(2.  2-10) 
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where  Kz  =  Kz(x-£a  •  y-na>  z- t,  ,w,  M)  and  Ky  =  Ky(x-£a  ,  y~*l  , 
s  S  ,i  rs  s* 

z-t,  ,w,  M)  are  evaluated  at  the  rcn  receiving  point  whose  coordinates  are 

cl 

given  by  x,  y,  z. 


Equation  (2.  2-8)  may  be  written  in  a  partitioned  form  as 


where  [  F  ] 


I ACpz ! 

(  AdPY  I 


(2.  2-11) 


The  boundary  conditions  at  the  lifting  surfaces  are  obtained  from 
the  motions  of  the  surfaces.  Specifically,  the  substantial  derivative  of  the 
surface  deflection  produces  the  required  normalwash  velocity  w  ,  Both 
lifting  surface  and  body  axis  singularities  are  introduced  to  produce  this 
normal  velocity.  If  lifting  surfaces  alone  are  considered  the  matrix 
formulation  is  taken  from  Equation  (2.  1-4),  i.  e.  , 


H  ■  w  {acp} 


When  axial  singularities  are  introduced  as  in  the  present  method  this 
formula  becomes; 


{»}  =  fDi{acP}  +  [F1  {acP(n} 


(2.  2-12) 


In  this  matrix  equation  only  the  boundary  conditions  on  the  lifting  surfaces 
are  considered.  The  matrix  [D]  represents  the  effect  of  lifting  surfaces  on 
themselves  and  is  square.  The  matrix  [F  ]  represents  the  effect  of  the  body 
axial  singularities  on  the  lifting  surfaces  and  is  rectangular. 
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Using  equation  (2.2-11)  in  (2.2-12)  gives 


(wrU 

:  <  W-W  \ 

= 

l  / 

l  B  / 

[D)  {*Cp} 


(2.  2-13) 


This  is  the  final  matrix  eqution  to  be  solved  for  the  unknown  lifting  pressures 

K )  • 

2.  2.  2  Interference 

The  slender  body  theory  satisfies  the  boundary  conditions  on  bodies 
undergoing  general  oscillatory  motions.  That  is,  the  slender  body  theory 
satisifies  the  upwash  and  sidewash  boundary  conditions  given  by  w  .  When 
lifting  surfaces  (or  other  bodies)  are  introduced  into  the  flow,  the  body 
boundary  conditions  are  no  longer  satisified.  The  region  of  greatest  dis¬ 
turbance  is  near  the  body-lifting  surface  intersection.  An  initial  attempt 
to  account  for  this  interference  is  presented  here  (Part  I).  The  basic  idea 
is  to  apply  lifting  surface  elements  directly  to  the  body  surface  in  regions 
where  the  body  boundary  condition  is  most  seriously  disturbed.  Each  lifting 
surface  element  possesses  its  own  wake.  Each  wake  trails  straight  back 
from  the  element  in  the  x-direction  (see  Reference  10).  Because  of  this 
fundamental  characteristic  the  configurations  that  may  be  considered  are 
limited.  Specifically,  placing  lifting  elements  on  the  actual  body  surface* 
is  not  acceptable  since  the  wakes  may  thread  in  and  out  of  the  body  surface 
aft  of  the  lifting  element  in  question  (see  Sketch  2.  2-1').  The  actual  wake,  if 
any  exists  at  all,  springs  from  the  separation  regions  of  the  body. 


*  Consideration  of  forward  or  rearward  inclined  surfaces  requires  a 
kernel  of  the  type  developed  by  Berman 
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Sketch  2.2-1 

The  shape  of  that  section  of  the  body  on  which  lifting  surface  elements 
are  placed  must  be  idealized.  Specifically,  the  interference  region  of  the 
body  must  have  a  constant  cross-sectional  shape.  Sketch  2.  2-2  presents  a 
graphical  example. 


Sketch  2.  2-2 


The  body,  then,  is  made  up  of  two  separate  components,  (1)  an 
axial  pressure  potential  doublet  distribution  whose  strength  is  obtained 
from  slender  body  theory  and  {2}  lifting  surface  elements  placed  on  the 
body  surface  whose  strength  is  obtained  by  a  solution  of  equation  (2,2-13). 


2.  3  The  Normalwash  Boundary  Conditions 

The  normalwash  boundary  conditions  for  the  lifting  surface  elements 
are  determined  by  assuming  the  surfaces  move  normal  to  themselves  with 
a  displacement  distribution  given  by  h.  These  normal  motions  are  given 
by  a  series  of  mode  shapes,  f.  and  generalized  coordinates  q. 

h  =  sX  fjqj  (2.3-1) 

The  normalwash  boundary  condition  is  given  by  the  substantial  derivative 
of  h,  i.  e.  ,  Dh/Dt.  The  normalwash  for  mode  i  is: 


Uco 


w. 

l 


d  (x/s) 


LJ  $ 


(2.  3-2) 


The  mode  shapes  may  be  input  numerically  (both  L  and  df./d(x/s)  may  be 
input)  or  as  polynomials. 


f. 

i 


5  X 
L  £ 

n=0  m=0 


a. 

mm 


(2.  3-3) 


Here  t  is  a  lateral  parameter  representing  the  radial  distance  from  the 
x-axis  to  a  point  on  the  surface  or  the  radial  distance  from  the  inboard 
edge  of  a  particular  lifting  surface  panel  to  a  point  on  the  surface  (see  Input- 
Procedure,  Part  I,  Vol  II  )•  The  constant  coef  is  simply  a  scale  factor. 
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The  x-derivative  of  f.  is  used  in  Equation  (2.  3-2)  and  is; 

5  5  /  \n“l  /  \m 

df/d(x/s)  =  coef  E  T, '  ainm  (2.3-4) 

n=0  m=0  's/ 

The  upwash  and  sidewash  to  the  bodies,  i.  e.  ,  are  also  determined  from 

Equation  (2.  3-2)  except  that  there  are  separate  modes  for  the  z  and  y-directions 
of  motion. 

The  slender  body  theory  of  Miles  requires  the  derivative  of 

dw  ^  _  _1_  dw 

dx  S  d(x/s) 

where 


d2(x/s) 

2  2 

If  modes  are  input  numerically  (f.  and  df./d(x/s)  then  d  f./d  (x/s)  must  be 

determined  numerically.  An  alternate  procedure  may  be  suggested  here. 

2  2 

Instead  of  f.  and  df./d(x/s)  the  terms  f.  and  d  f./d(x/s)  could  be  input.  The 
i  i  2  2 

term  df./d(x/s)  is  then  obtained  by  numerical  integration  of  d  R/d(x/s) 

which  is  more  accurate  than  numerical  differentiation*. 

Also  required  for  the  slender  body  solution  is  dRQ/dx.  This  quantity 
may  be  obtained  numerically  from  the  input  Rq.  An  alternate  approach  may 
again  be  suggested  whereby  dRQ/dx  is  input  and  Rq  is  obtained  by  numerical 
integration. 

*  Actually  the  first  value  of  dT/d(x/s)  is  input  along  with  the  second  derivatives 
to  establish  the  proper  level  for  df./d(x/s). 


1  I  d“fi 


d2(x/s) 


<1 


i«.\ 


df. 


I  d(x/s) 


(2.3-5) 


5  5 


E  E  (I)m 

n=0  m=0  '  s  '  v  s  ' 
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2.4  Loads  and  Forces 


2.  4.  1  Generalized  Forces 


The  lifting  pressures,  AC  ,  on  the  lifting  surfaces  are  given  by  the 

P  m 

solution  of  (2.2-13).  The  lifting  pressures,  AC  1  on  the  bodies  are  given 

P 

by  (2.2-4).  These  two  distributions  may  be  combined  into  one  partitioned 

array,  AC  as  follows: 

P 


(  ACd  ) 

1  AC  1 

— 

K,ty! 

1  p  1 

(2.4-1) 


The  lifting  pressure  AC  is  composed  of  two  components,  i.  e.  ,  AC„ 

P  A 

and  ACy. 


AC 


(f) 


AC, 

l 

AC, 


The  lifting  pressure  AC^  acts  on  the  area  of  the  body  projected  onto  the 
z  =  0  plane.  The  lifting  pressure  ACy  acts  on  the  area  of  the  body  projected 
onto  the  y  =  0  plane.  All  loadings  are  in  the  form  of  pressures  acting  on 
areas  and  the  calculation  of  the  generalized  forces  may  proceed  on  this  basis. 
There  are  as  many  lifting  pressure  distributions  as  there  are  modes.  The 
subscript  "i"  identifying  the  mode  will  be  added  to  AC^. 


AC 


2  ACp  q. 


The  generalized  forces,  Q„,  associated  with  the  generalized  coordi¬ 


nates,  q.,  are  determined  using  the  principle  of  virtual  work, 

^  fch 

work,  6Wj  associated  with  the  j  mode  is: 


The  virtual 


6W. 


//  * 


AC 


6h  ds 


(2.4-2) 


L.  S. 


The  sane  virtual  work  is  obtained  using  the  generalized  forces 


6  W, 
) 


-2  s' 


E 


% 


6q. 


(2  4-3) 


Recalling,  from  (2,3-1),  tl  at 


and  equating  (2.  4-2)  and  (2,  4-3)  gives 


-2 


2 

s  q 


If 


i. 


ds 


Solving  for  CE.  gives: 


(2.  4-4) 


Numerically,  the  double  integral  is  replaced  by  a  summation: 


where  As  is  the  box  area  on  the  lifting  surfaces  and  where  As  is 
Rq  A£n  on  the  projected  body  surfaces.  The  factor  1/2  is  accounted  for 
by  restricting  the  range  of  the  subscript  "n"  to  the  right-half  of  the  aircraft. 
(Even  the  projected  area  of  the  body  is  halved,  i.  e.  ,  As^  =  RQA£;n  if  it 
lies  on  the  plane  of  symmetry.  ) 

2.4.2  Aerodynamic  Parameters 

It  is  desirable  and  sometimes  necessary  to  generate  conventional 
aerodynamic  data.  Such  data,  in  addition  to  being  useful  in  Lbemselves, 
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provide  an  excellent  check  for  the  computer  program  and/or  specific  cases 
to  be  run  by  it.  For  example,  a  check  of  the  span  loading  for  rigid  body 
pitching  and/or  plunging  could  bring  to  light  some  geometrical  data  error 
in  a  particular  case.  Thus  the  aerodynamic  data  output  acts  as  a  program 
monitor. 

The  local  normal  force  coefficient  and  pitching  moment  coefficient 
about  the  local  leading  edge  is 


*/ 


AC  d£ 
P 


(2.  4-6) 


c 2  l 


AC  (I  -  x  )  dg 


(Z  4-7) 


where  c  is  the  local  chord  length.  The  total  lift  and  side  force  c(  efficients 


(1  +6) 


c  c  d1)  + 
n 


J  RoACzdl 


(2.4-8) 


(1  -6) 


c  c  d£  +  g 
n  b 


/  Ro 


ACy  d£ 


(2.4-9) 


Here  A  is  the  reference  area,  6  is  a  symmetry  flag  (6=1,  for  symmetry, 
6=  -1,  for  antisymmetry  and  6=0  for  asymmetry).  The  limit,  R.  S. ,  indicates 
that  the  integration  is  to  be  carried  out  over  the  right  side  of  the  aircraft  and 
the  limit  B.  indicates  that  the  integrals  are  to  be  taken  along  the  lengths  of  the 
bodies.  The  term  g  is  defined  in  Eq.  (2.  4-12).  The  pitching  moment,  yawing 
moment  and  roiling  moment  coefficients  are  given  by; 


(1  +6) 


( c  c  -  (xT  _  -  xn)cc  )  dO 

m  JL.E.  R  n 


(2.  4-10) 


-g  J  ACz({-xR)Rod«  (  +  nose  up) 
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where  x^  is  the  point  about  which  moments  are  taken.  The  rolling  moment 
is  taken  about  the  x-axis. 

Dynamic  stability  derivatives  can  be  obtained  from  the  complex 
coefficients  just  described.  Each  complex  aerodynamic  coefficient 
possesses  an  in-phase  and  an  out-of-phase  component,  the  real  part  being 
the  component  in  phase  with  the  motion  and  the  imaginary  part  the  component 
out  of  phase  with  the  motion  (90°  out  of  phase).  Reference  11  describes  a 


26 


method  for  determining  dynamic  stability  derivatives.  Specifically,  if  C 
represents  any  aerodynamic  coefficient  such  as  cn  ,  c  ,  C  ™  C  ^  etc.  , 
and  if  the  subscripts  a  and  6  represent  modified  plunging  and  pitching  then:* 


C  a  ~ 

C* 

(kr  =  0) 

static  derivative 

Oa  = 

1 

k 

r 

Im  Cq,  (kr  =  «) 

damping  derivative 

c« 

1 

',2 

Re  {  Cq  (kr  =  0)  - 

Cq  (k^  =  e)}  acceleration  derivative 

The  program  must  be  run  at  two  different  values  of  reduced  frequency 
=  0  and  k^  =  «  to  obtain  these  dynamic  derivatives.  The  value  of  < 
is  usually  small  (about  0.  1).  The  dynamic  derivatives  will  then  be  valid 
for  aircraft  motions  in  the  frequency  range  0  to  € . 

The  values  of  the  q  or  pitch-rate  derivatives  are  obtained  as  follows; 

1 

Cq  =  ~  Im  ^6  ^r  "  ^ 

M  r 

1 

C-  =  —  Re  {  C0{ kr  =  0)  -  Co  (kr  =  c  )}-  C-a 

q  [<- 


Similar  expressions  for  Cp,  Cp,  Cp,  and  are  obtained  if  the 
values  of  Cq  and  Cq  are  replaced  by  values  for  modified  sideslipping  and 
yawing.  The  expressions  for  the  rolling  damping  and  acceleration 
derivatives  are  obtained  as  follows: 


where  C^  is  the  complex  aerodynamic  coefficient  due  to  unit  rolling  velocity 

*  The  term  "modified"  indicates  that  the  amplitude  of  translation  (plunging 
or  sideslipping  is  normalized  by  the  term  (ik^).  The  resulting  derivatives 
then,  pertain  to  a  and  P  respectively. 


27 


2.  5  Aerodynamic  Influence  Coefficients  and  Harmonic  Gust  Coefficients 


2.5.1  General  Considerations 

The  aerodynamic  influence  coefficient  matrix  used  in  this  method  is 
described  in  References  12  and  13.  The  basic  idea  is  to  produce  a  set  of 
influence  coefficients  that  are  independent  of  mode  shape  and  smaller  in 
number  than  the  number  of  elements  or  boxes.  An  obvious  set  of  AIC's 
would  be  [  D  ]  *  whereby  lifting  pressures  are  obtained  directly  from  the 
normalwash  boundary  conditions.  A  second  example  would  be  the  matrix 

of  generalized  forces  f  Q],  a  typical  element  of  which  is  Q...  The  advantage 

_  1  H 

of  [  D  ]  is  that  it  is  independent  of  mode  shape;  the  advantage  of  [Q]  is  that 

it  is  small  and  easily  handled.  A  compromise  between  f  D]  *  and  [Q]  is 

the  current  [  AIC  ].  A  preselected  set  of  "submodes"  is  selected  for  each 

lifting  surface  strip  or  bay.  The  submode  allows  motion  only  on  the  strip 

or  bay  considered.  A  set  of  generalized  forces  is  obtained  from  each  of 

these  submodes.  For  instance  if  the  submodes  are  pitching  and  plunging 

and  there  are  15  bays  then  the  AIC  will  be  a  30  x  30  matrix  of  generalized 

forces.  It  is  assumed  that  the  general  mode  shapes  may  be  built  up  of 

these  submodes  by  superposition.  The  generalized  force  matrix  for  any 

specific  mode  shape  is  then 

Qjj  =  { )T  [AIC]  {  f.  } 

where  f.  are  components  of  the  i^  mode  shape  in  terms  of  the  submodes. 

The  value  of  using  AIC's  lies  in  the  fact  that  once  the  matrices  of 
AIC's  have  been  generated  for  a  specific  planform,  Mach  number  and 
reduced  frequency,  any  number  of  aeroelastic  analyses  can  be  performed 
in  which  only  the  inertial  and  stiffness  properties  are  altered,  since  the 
AIC's  are  independent  of  the  aircraft  vibration  mode  shapes  and/or  static 
deflection  modes. 
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The  aerodynamic  influence  coefficients,  as  derived  in  Reference  1?., 
relate  the  oscillatory  aerodynamic  moments  and/or  forces  acting  at  the 
specified  AIC  control  points  to  the  harmonic  rotations  and/or  deflections 
of  these  control  points.  By  the  definition  of  Reference  12  the  equation  for 
the  dimensional  AIC  matrix  is 

i  U°°  i 

tchld  =  2  Po  ~  {ZI2)  fBl  ED1  fw1  (Z  5"1} 

r 

Uco 

where c/2is  the  reference  semi -chord,  ~r -  is  the  reference  reduced 

br  u> 

velocity,  [  B  ]  is  the  integration  matrix  described  in  Sec.  2.  5-3,  and 
[W  ]  is  the  substantial  derivative  matrix  described  in  Sec.  2.  5-2. 

A  set  of  dimensional  static  AIC's  [C^s^d  rnay  also  be  defined  by 

[Chs]d  "  '2  P0(?/2)[B]  [D  (kr=0)]-’  [W]  (2.5-2) 

These  can  be  used  in  static  aeroelastic  analyses  of  lift  effectiveness, 
divergence,  control  surface  effectiveness,  reversal,  and  so  forth. 

The  harmonic  gust  coefficients  were  also  derived  in  Referenced 
for  a  gust  in  the  plane  of  the  surface.  For  a  gust  field  in  a  plane  with 
arbitrary  dihedral*  rg,  Eq.  (39)  of  Reference  12  is  generalized  to  read 

{c}  =  j  po  U*  wg  [  B  ]  [  D  ]" 1  {  cos  (rg-Y)  exp  (-ikr  ^^-)}(2.5-3) 

where  U<x>  is  the  aircraft  velocity,  W  is  the  harmonic  gust  amplitude, 

O 

is  the  reduced  frequency,  Y  is  the  local  surface  dihedral,  and  xq  is 
the  gust  reference  coordinate. 


*  Vertical  velocities  are  obtained  when  the  dihedral  angle  is  zero.  If 

the  dihedral  angle  is  -90°  the  gust  velocities  are  in  the  lateral  (+y)  direction. 
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2.  5.  2  The  Substantial  Derivative  Matrix 

The  Doublet-Lattice  Method  requires  the  determination  of  the 
oscillatory  normalwash  on  the  lifting  surface  at  the  3/4-chord  of  each  box 
In  the  derivation  of  the  substantial  derivative  matrix,  [W  ],  we  assume  that 
the  aerodynamic  idealization  for  the  particular  aeroelastic  problem  pro¬ 
vides  a  sufficient  number  of  (spanwise)  strips  on  the  lifting  surface  so  that 
the  oscillatory  normalwash  on  each  box  can  be  determined  only  by  chord- 
wise  interpolation  and  differentiation  of  the  assumed  local  deflection  mode 
shapes  corresponding  to  each  local  aerodynamic  degree  of  freedom.  It  is 
assumed  that  an  interpolation  matrix,  [Wj],  can  be  found  such  that 

[l>3 C/41  "  Iwll  t  h  1 

where  I13 c/4  denotes  the  deflection  at  the  3 /4-chord  point  of  each  box  and 
h  denotes  the  AIC  control  point  deflections.  The  matrix  [Wj]  has  the 
partitioned  form 


(2.  5-1) 


where  the  superscripts  indicate  strips  in  the  aerodynamic  idealization  of 
the  lifting  surface.  The  oscillatory  normalwash  at  the  3/4-chord  of  each 
box  is  given  by 


w(x) 

u® 


dh(x) 

dx 


h  (x) 


(2.  5-5) 


Generalizing  Eq.  (2.  5-5)  in  matrix  notation  we  have 


where 


{w}  =  [W  ]  {h} 


[Wl 


dW 

dx 


k 

+  5  — —  [W  1 

c/2  1 2 


(2.  5-6) 


(2.  5-7) 


Several  alternatives  to  the  computation  of  the  [W  ]  matrix  can  be 
used  depending  on  the  manner  in  which  the  deflection  characteristics  of 
the  aircraft  are  described.  Ref.  12  describes  three  possible  alternatives 
for  the  [W  ]  matrix;  Alternative  No.  1  is  used  in  the  generation  of  the  AIC's 
for  planar  surfaces  by  computer  program  H7WA.  The  nonplanar  computer 
programs,  H7VVB  and  H7WC  have  the  capability  of  generating  AIC's  for  non¬ 
planar  surfaces  according  to  either  of  the  two  alternatives  of  the  [W  ]  matrix 
described  in  Reference  13  in  Appendices  I  and  II,  which  are  reproduced  in 
Appendix  D. 

The  two  alternatives  describe  two  types  of  submodes  that  can  be  used. 
The  submodes  for  each  bay  are: 

1)  Alternative  ft 1 :  plunging,  pitching,  control  surface  rotation  and 
tab  rotation  or  control  surface  plunging. 

2)  Alternative  #2.  Three  cambering  modes,  control  surface  rotation, 
or  tab  rotation. 

Alternative  #1  makes  available  six  degrees  of  freedom  or  submodes 
per  bay,  while  Alternative  §2  furnishes  seven.  Of  these  available  degrees 
of  freedom  or  submodes  only  four  may  be  used  at  one  time.  The  reason 


for  this  restriction  is  the  limited  core  capacity  of  the  computer.  More  of 
these  available  submodes  may  be  made  available  if  a  core  larger  than  32K 
is  available. 


where  the  h^4  are  the  deflections  at  the  1 /4-chord  of  each  box  and  h  are 
again  the  AIC  control  point  deflections,  f  Z  ]  has  the  same  partitioned  form 
as  the  matrix  [  W ^  ]  (Eq.  2.  5-4). 
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The  generalized  forces  acting  at  the  1 /4-chord  ol'  each  box  are 
given  by 


{  Fc/4  }  =  q  |" A J  { A C p }  (Z.  5-10, 

where  |\aJ  is  the  diagonal  matrix  of  box  areas.  By  applying  the  principle 
of  virtual  work  it  is  seen  that 

(F  }AK  =  [Z]T  {Fc/4}  (2.5-M) 

By  combining  Eqs.  (2.  5-1  0)and(2.  5-1 1)  and  comparing  with  Eq.  (2.  5-8)we 
have 


[B]  =  [Z]T  [aJ  (2  5-12) 

The  actual  form  of  the  [B  ]  matrix  depends  on  the  manner  in  which  the 
deflection  characteristics  of  the  aircraft  are  described,  and  hence,  on  the 
form  of  the  [  W  ]  matrix.  The  Lwo  alternatives  for  the  [  B  1  matrix  are 
described  in  Reference  12  along  with  the  corresponding  [W  ]  matrix 


2.  6  Matrix  Solution 


When  modes  are  considered  the  matrix  equation  (2.  2-13)  must  be 
solved  for  each  mode.  The  symbolic  solution  for  multiple  modes  is 


[  A  Cp  ]  -•  [  D  ]"  [  w  *  wB  ] 


(2. 6-1) 


where  the  matrices  [  w  -  ]  and  [ACp  )  are  the  normalwash  and  resulting 

pressures  for  all  of  the  modes.  Each  column  of  either  [A  Cp  1  or  [  w  -  ] 

represents  one  mode.  For  AIC's  Equation  (2.  5-1)  must  be  solved- 


tCiJd  =  2  Po  “lT-(c/2)fB]  fD)"1  [W] 

r 


(2  6-2) 


In  both  cases  the  symbolic  solution  contains  a  term  like  f  D  J  f  wr  ]  where 
[  wr  ]  is  [  w  -  Wg  ]  in(2.  6-l)and  [  W  ]  in  (2.  6-2).  Define  the  term  f  D )  ^ 

[  wr  ]  as  [  P  ]. 


[P]  =  [Df1  [wr] 


(2.6-3) 


In  the  present  method  the  desired  solution  matrix  [P]  is  not  obtained 
exactly  as  shown  in  this  equation.  The  inverse,  in  the  usua-1  sense  of  the 
word,  is  not  determined  since  this  would  be  very  inefficient  for  large  matrices 

The  solution  matrix  [P]  is  obtained  directly  from  the  equation 

[  wr  ]  =  [  D  ]  [  P  ]  (2.  6-4) 


by  direct  Gaussian  Triangularization  and  back  solution.  This  subroutine  is 
called  SOLVIT  and  is  described  in  Reference  14,  The  disadvantage  of  this 
procedure  is  that  future  solutions  can  be  obtained  only  by  a  repetition  of  the 
triangulariv.ation  and  back  solution  (a  (1/3)N  +  N  M  process),  where  N  is 

the  order  c„  the  matrix  [  D  ]  and  M  is  the  order  of  [wr  ]  ).  Also  the  solution 
matrix  [  P  ]  must  fit  into  the  computer  core.  When  [  P  ]  is  small,  as  with  the 


34 


modal  approach,  this  restriction  is  not  important.  However  for  AIC's 
the  solution  matrix  may  be  larger  than  the  available  storage. 

A  generalization  of  this  method  may  be  made  to  eliminate  both 

disadvantages  discussed  above  (see  Reference  14,  page  37).  The  basic 

idea  is  to  save  the  instructions  needed  to  triangularize  the  matrix  [D]. 

It  just  so  happens  that  these  instructions  come  in  the  form  of  constants,  the 

number  of  which  fit  into  the  empty  space  generated  in  [  D  ].  by  triangulari- 

zation.  The  triangularized  form  of  [D]  fits  into  the  upper  triangular  sector, 

including  the  diagonal,  while  the  constants  of  triangularization  fit  into  the 

lower  triangular  sector,  excluding  the  diagonal.  This  matrix  is  termed 

the  QuasiTnverse  since  both  halves  of  the  matrix  can  be  used  to  find  a 

2 

solution  using  only  N  operations  per  mode  or  submode.  The  Quasi-Inverse 
requires  identically  the  same  number  of  operations  as  a  direct  solution, 

3 

i.  e.  ,  (1  / 3) N  operations.  In  fact  this  modified  form  of  SOLVIT  works 
exactly  like  SOLVIT  if  the  solution  matrix  fits  in  core  (except  that  the 
Quasi-Inverse  matrix  must  be  read  on  tape  as  it  is  generated  so  that  it  can 
be  used  in  the  future.  ) 

The  difficulty  of  running  future  solutions  is  thus  solved  Future 

solutions  are  solved  simply  by  calling  back  the  Quasi-Inverse  and  performing 
2 

an  N  M  operation  (where  M  is  the  number  of  modes  or  submodes). 

The  restrictive  requirement  that  the  solution  matrix  fit  into  core  is 

also  eliminated  since  the  normalwash  matrix  may  be  split  into,  say,  two  or 

more  parts,  each  of  which  fit  into  the  computer  core.  The  first  part  of  the 

solution  corresponding  to  the  first  part  of  the  normalized  matrix  obtained 

3  2 

during  the  first  pass  through  the  computer((l /3)N  +  N  M  operations). 

Also  formed  during  this  pass  is  the  Quasi -Inverse.  The  second  and  sub¬ 
sequent  passes  through  the  computer  generate  the  remainder  of  the  solution 
2 

(MN  operations  each).  The  second  and  remaining  parts  are  thus  considered 
future  solutions.  This  method  is  programmed  in  two  parts  QUAS  and  FUTSOL 
QUAS  forms  the  Quasi-Inverse  and  solves  for  the  first  part  of  the  solution 
and  FUTSOL  generates  future  solutions  from  the  Quasi-Inverse. 
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3.  0  CALCULATED  RESULTS 


3.  1  Convergence  of  Results 

The  number  of  boxes  necessary  to  give  a  converged  result  is 
dependent  on  the  reduced  frequency.  For  the  steady  case  the  convergence 
is  very  rapid.  In  fact,  one  box  chordwise  gives  usable  results  (the  method 
of  Weissinger  uses  one  box  per  spanwise  station). 


For  the  other  extreme,  many  boxes  are  necessary  for  very 
high  reduced  frequencies.  As  an  example, Figure  1  presents  a  convergence 

m2  g 

study  for  a  reduced  frequency  (based  on  the  span  2s,  k  =  — )  of  1. 5. 

r  l  Vqo 

As  shown  in  the  figure  the  wave  length  corresponding  to  this  frequency  is 
approximately  the  length  of  the  wing-tail  assembly.  This  means  that  the 
effect  of  the  wing  on  the  tail  is  out  of  phase  with  the  wing  motion  about  '80°. 
Thus  the  downwash  of  the  wing  on  the  tail  is  opposite  to  the  wing  motion. 
This  is  a  very  severe  case  and  many  chordwise  boxes  are  required.  The 
figure  shows  that  approximately  25  boxes  per  wave  length  are  required. 
Thus  for  this  case  the  box  length  AXgis  related  to  the  reduced  frequency  by: 


Ags.  < 


13k. where  kr 


wc 


Actually  the  wing  and  tail  box  lengths  are  not  equal  for  Figure  1. 
The  number  of  chordwise  boxes  for  wing  and  tail  are: 


TOTAL 

WING 

TAIL 

9 

5 

4 

11 

6 

5 

14 

8 

6 

18 

10 

8 

22 

12 

10 

There  are  twelve  spanwise  strips  for  the  case  presented  in  Figure  1.  A 
case  using  eight  spanwise  strips  is  also  shown  for  the  case  of  22  total 
chordwise  elements.  The  error  for  the  real  and  imaginary  parts  is  about 
2.  5%  and  0.  80%,  respectively. 


3.2  Nearly  Coplanar  Wing-Tail  Configurations 


The  difficulty  encountered  by  the  Doublet-Lattice  Method  of 
Reference  1  for  nearly  coplanar  wing-tail  combinations  is  discussed  in 
Section  2.  1.  The  refinement  of  the  method  is  also  presented  there.  Figure  2 
presents  results  for  a  nearly  coplanar  wing-tail  configuration  using  both  the 
original  (unrefined)  and  the  refined  methods. 

The  wing-tail  assembly  is  made  to  plunge  with  a  reduced 
frequency  of  0.  6  (based  on  a  chord  length  of  0.  87424  s).  The  total  lift 
coefficient  (wing  plus  tail)  is  given  (based  on  total  area  of  both  surfaces). 

It  can  be  seen  that  the  original  method  starts  to  lose  accuracy  when  the  gap 
to  chord  ratio  is  about  0.  1.  (Notice  that  the  widths  of  the  first  eight  strips 
is  0.  1).  The  refined  Doublet-Lattice  Method  gives  accurate  results  for 
small  z/c.  Appendix  B  explain0  that  a  numerical  limit  must  be  taken  in 
the  new  or  refined  method.  Any  inaccuracies  generated  from  this  limit  are 
not  apparent  even  at  a  value  of  z/c  of  0.  025.  Thus  the  refined  method  is 
valid  for  all  practical  cases  of  interest. 

The  wing-tail  configuration  of  Figure  1  is  used  to  present  a 
.-.econd  "gap"  study  (see  Figure  3).  A  much  wider  range  of  vertical  separa¬ 
tion  (gap)  between  the  wing  and  tail  is  used?  Also  the  reduced  frequency  is 
considerably  higher.  Again  the  assembly  is  made  to  plunge  and  the  total  lift 
coefficient  (based  on  the  area  of  the  wing  alone)  is  given  versus  vertical 
spacing  or  "gap".  Although  not  entirely  evident  the  curves  of  real  and 
imaginary  parts  of  the  lift  coefficient  have  horizontal  tangents  as  z— 0  since 
they  are  symmetrical  with  z. 

3.  3  Wing-Body  Combinations 

The  methods  of  Woodward  (Ref.  6)  and  Giesing  (Ref.  3)  handle 
wing-body  combinations  in  steady  flow.  A  recent  unpublished  extension  of 
Reference  3  allows  bodies  of  varying  diameter  to  be  handled.  These  methods 
can  be  used  to  furnish  test  cases  for  the  present  method,  at  least  for  the 
static  case. 


#  The  tail  moves  aft  as  it  moves  up  along  the  line  x  =  0,  875  z  . 


Figure  4  presents  the  span  loading  (cc^/c)  for  a  wing -fuselage 
combination  as  calculated  by  the  present  method  and  the  methods  of 
Woodward  (Ref.  6)  and  Giesing  (modified  Ref.  3).  The  configuration  shown 
in  the  figure  consists  of  a  low  aspect  ratio  wing  attached  at  midplane  to  a 
thick  parabolic  body  of  revolution.  The  center  of  the  wing  root  io  at  the 
center  of  the  fuselage.  The  maximum  diameter  is  1/3  the  span.  The  con¬ 
stant  section,  upon  which  lifting  surface  elements  are  placed,  is  half  the 
length  of  the  body  and  centered  longitudinally.  The  span  loadings  shown  in 
Figure  4  are  entirely  fuselage  induced  since  the  wing  is  held  at  zero  incidence 
In  Figure  4(a)  the  fuselage  is  given  an  angle  of  attack  of  1 . 0  radian.  It  can 
be  seen  that  the  present  method  agrees  almost  exactly  with  the  method  of 
Woodward*  since  the  idealization  for  these  two  methods  is  about  the  same. 
The  disagreement  with  the  method  of  images  (modified  version  of  Ref  3) 
is  due  to  the  fact  that  a  different  idealization  is  used.  The  boundary  condi¬ 
tion  on  the  fuselage  surface  is  satisfied  using  a  system  of  images,  one  for 
each  horseshoe  vortex,  plus  an  axial  singularity  system. 

The  span  loading  across  the  fuselage  is  induced  by  the  wing 
loading  (which  was  induced  by  the  body  incidence  or  camber)  This  loading 
is  the  fuselage  lift  carry-through.  Slender  body  theory  predicts  no  net 
load  since  the  body  closes. 

Figure  4(b)  shows  the  span  load  on  the  wing  whe„n  the  fuselage 
possesses  a  parabolic  camber  with  a  maximum  deflection  equal  to  the 
maximum  radius.  The  agreement  between  the  present  method  and  Woodward' 
method  is  almost  exact. 

Figure  5  presents  a  comparison  of  calculated  spanwise  load 
distribution  obtained  using  three  slightly  different  approaches  to  the 
interference  problem.  The  configuration  and  flow  condition  used  is  that  of 
Figure  4(a).  The  circles  represent  the  present  method  and  are  the  same  as 
those  in  Figure  4(a).  The  squares  represent  a  method  that  differs  onlv 
slightly  from  the  present  method.  The  subtle  difference  between  the  two 
methods  is  that  the  flow  field  due  to  '  body  axial  singularities  is  calculated 


*  The  control  point  must  be  placed  at  85%  of  the  box  to  render  Woodward's 
results  correct.  See  Reference  5. 


on  all  lifting  surface  panels,  including  the  body  panels,  and  not  just  on  the 
wing  panel  as  in  the  present  method.  Since  the  body  panels  do  not  lie  exactly 
on  the  body  surface, the  total  normalwash  there  is  different  from  zero  (since 
the  axial  singularities  were  meant  to  cancel  the  normal  flow  only  at  the 
exact  location  of  the  body  surface).  This  slight  difference  is  enough  to 
cause  large  differences  in  the  spanwise  loading  on  the  body  surface  even 
though  the  loading  on  the  wing  compares  favorably  with  the  present  method. 

A  second  alternative  (method  of  Reference  5)  is  to  use  body  panels 
alone  to  simulate  the  body  surface.  For  this  method  no  axial  singularities 
are  used.  When  the  body  possesses  no  angle  of  attack  (boundary  condition 
of  zero  normalwash  to  the  body  surface)  the  results  are  accurate  (see 
Reference  5,  Figure  12).  However, when  the  body  possesses  an  angle  of 
attack  or  camber, then  the  results  seem  to  be  very  inaccurate  especially 
over  the  body  surface.  Figure  5  shows  such  a  calculation.  The  reason  for 
this  inaccuracy  at  angle  of  attack  is  that  the  body  panels  are  simulating  an 
annular  or  ring  wing.  Such  an  annular  wing  is  much  more  effective  in  lift 
than  the  body  it  represents.  If  the  body  has  no  region  of  separation  and  if 
no  vortices  are  shed  from  the  body  then  the  only  lift  it  can  carry  is  the  lift 
carry-through  which  is  induced  by  the  wings.  The  spanwise  distribution  ol 
lift  carry-through  has  the  characteristic  shape  shown  on  Figures  4(a'i  and  4(b) 
The  lift,  over  and  above  the  lift  carry-through,  is  caused  by  the  annular 
wing  effect  of  the  body  elements.  A  possible  solution  to  this  problem 
would  be  to  use  wakeless  elements  on  the  body  surface.  Such  wakeless 
elements  can  be  generated  using  two  steady  or  unsteady  horseshoe  vortices: 
one  placed  behind  the  other  such  that  the  wakes  cancel  out. 

Figure  6  presents  an  example  of  an  unsteady  wing-fuselage 
calculation.  The  configuration  used  is  that  of  Figure  4.  The  two  modes 
of  oscillation  are  again  the  same  as  those  used  in  Figure  4,  for  the  static 
case.  Specifically: 

Mode  1  (Figure  6(a))  z  =  (x  -  xc)  eltvf 

x  -  x  , 

Mode  2  (Figure  6(b))  z  =  4.  0  R  ( —  — )  elw 

ITJcLX  T  T 
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Here  x  ,  L  and  R  are  the  center,  length  and  maximum  radius  of 
c  max 

the  body  respectively.  The  frequency  of  oscillation  (based  on  the  average 
chord,  c"=  7)  is  =  1.0.  The  wing  remains  stationary  while  only 

the  body  is  in  motion.  The  span  loading  across  the  wing  is  induced  by  the 
body  motions  while  the  span  loading  across  the  body  is  induced  by  the  loading 
on  the  wing.  The  span  loading  across  the  fuselage  does  not  include  the 
unsteady  slender  body  component*  and  is  just  the  lift  carry-through. 

A  second  wing-body  case  is  presented  in  Figures  7,  8,  and  9.  The 
specific  configuration  considered,  which  is  shown  m  Figure  7,  is  a  wing 
with  an  engine  nacelle  mounted  at  its  tip.  The  wing  nacelle  combination  is 
made  to  plunge  and  pitch  harmonically  with  a  reduced  frequency,  (based  on 
the  semispan),  of  1.72.  The  pitch  axis  is  shown  in  Figure  7  The  real  and 
imaginary  parts  of  the  distribution  of  lift  coefficient  are  given  for  both 
pitching  and  plunging.  In  addition,  the  distribution  of  lift  coefficient  for  the 
wing  alone  is  presented  for  reference.  Notice  that  the  nacelle  has  a  sub¬ 
stantial  effect  on  the  distribution  of  lift  coefficient.  This  effect  arises  from 
two  sources;  the  first  is  the  upw'ash  generated  by  the  motion  of  the  nacelle 
while  the  second  is  the  end-plate  effect  of  the  nacelle  or  the  wing  tip. 

The  slender  body  simulation  of  the  nacelle  is  accomplished  using  two 
bodies.  The  axes  of  the  bodies  are  centered  within  the  nacelle  spanwise 
and  given  a  vertical  displacement  of  ±0.  108252.  The  radius  of  these  bodies 
is  0.  108252.  Even  though  these  two  bodies  do  not  quite  fill  up  the  nacelle 
cross-sectional  area,  it  can  be  shown  that  this  choice  of  center  and  radius 
gives  the  best  representation  of  the  flow  around  the  nacelle  cross-section 
(in  two-dimensions). 

Figures  8  and  9  present  comparisons  of  experimental  and  calculated 
lifting  pressure  coefficient.  Specifically,  Figure  8  gives  the  pressure  at 
two  spanwise  stations,  y/s  -  0,27  and  0.715  for  the  plunging  case.  The  wing- 


*  Unsteady  slender  body  theory  gives  a  non-zero  lifting  force  on  a  body 
even  if  a  closed  body  is  considered. 
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alone  pressure  distribution  is  also  given  at  the  spanwise  station  closest  to 
the  nacelle;  i.  e.  ,  y/s  =  0.  715. 

Figure  9  gives  the  pressures  at  the  same  two  stations  for  the 
pitching  case.  The  test  and  theory  correlate  best  at  station  y/s  =  0.  27  for 
both  the  pitch  and  plunge  case.  This  may  be  due  to  flow  separation  on  the 
aft  portion  of  the  wing  which  seems  to  be  present  in  the  vicinity  of  y/s  0.  716. 
The  loss  <_>f  lift  aft  of  about  the  80^  chord  point  indicates  that  this  is  the  case. 

3. 4  Conclusions 

(1)  The  Doublet- Lattice  Method  is  a  simple,  versatile  and 
accurate  lifting  surface  theory.  The  method  is  capable 
of  analyzing  lifting  surfaces  with  arbitrary  planform  and 
dihedral.  Control  surfaces,  cither  full-  or  partial-span, 
may  be  included.  Problems  of  intersecting  and/or  inter¬ 
fering  nonplanar  configurations,  such  as  a  wing-pylon 
combination,  a  T-  or  V-tail,  a  wing-tail  combination,  etc.  , 
may  be  analyzed.  The  method  documented  here  is  also 
capable  of  solving  problems  involving  lifting  surfaces  and 
bodies  where  the  bodies  may  be  in  motion  Included  in  the 
calculations  are  options  for 

(a)  Aerodynamic  data  including  lifting  pressures,  span- 
wise  lift  and  moment  distributions,  aerodynamic 
center  locations,  total  lift  and  side  force  coefficients, 
and  total  pitching,  yawing  and  rolling  moments; 

(b)  Generalized  forces  for  polynomial  mooes  of  motion 
specified  by  the  user; 

(c)  Aerodynamic  Influence  Coefficients; 

(u)  Gust  loads  from  '  harmonic  gust  field; 

(e)  Symmetry  and  ground  effects. 
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{2}  A  refined  version  of  the  original  Doublet-Lattice  Method 
is  required  for  nearly  coplanar  wing-tail  combinations 
The  refinement  consists  of  considering  the  planar  and  non- 
planar  parts  of  the  Kernel  separately.  The  details  are 
given  in  Appendix  B  As  a  control  point  approaches  the 
plane  of  the  sending  element  a  limit  is  taken  numerically 
This  limiting  procedure  is  accurate  for  all  practical  values 
of  vertical  displacement. 

(3)  The  number  of  chordwise  boxes  must  increase  in  proportion 

to  the  frequency;  specifically,  Axg/c  s  ir/25kr  Other  require¬ 
ments  exist  for  the  placement  of  strips  and  boxes  These 
requirements  are  discussed  in  Part  I,  Vol  II. 

(4)  The  boundary  condition  on  the  lifting  surface  elements, 
placed  on  body  surfaces  for  wing-body  interference  purposes, 
must  be  zero  normalwash.  The  elements  on  the  body  surface 
represent  a  ring  or  annular  wing*.  This  annular  wing  must 
be  used  for  interference  purposes  only;  any  angle  of  attack 
will  produce  an  undesirable  annular  wing  lift.  A  slender  body 
theory  for  circular  bodies  is  used  to  reduce  the  normalwash, 
generated  by  a  body  in  motion,  to  zero. 

(5)  The  tube  on  which  lifting  surface  elements  are  placed,  to 
generate  wing-body  interference,  must  have  a  constant  cross- 
section  along  its  length.  The  cross-sectional  shape  may  be 
noncircular.  The  reason  for  this  is  that  each  element  has  a 
wake  trailing  straight  back  downstream  to  infinity  If  the 
tube  were  of  varying  radius,  element  wakes  would  be  thread¬ 
ing  in  and  out  of  the  body  surface.  Such  an  idealization  has 
not  been  tried  and  its  validity  has  not  been  proven. 

The  synthesis  of  a  noncircular  slender  body  theory  out 
of  one  for  circular  bodies  is  accomplished  by  placing  more 
than  one  circular  body  in  close  proximity.  The  distance 

*  Actually,  the  body  cross-section  may  be  of  arbitrary  shape  even  though  the 
slender  body  theory  is  for  a  circular  cross-section.  A  slender  body  theory  for 
cross-sections  of  different  shape  could  replace  this  simple  theory. 
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between  axes  and  the  radius  of  the  various  bodies  may 
be  obtained  fron  a  two-dimensional  analysis.  This 
approach,  although  only  approximate,  is  reasonably 
accurate.  In  general,  the  interference  generated  by 
the  body  elements  is  less  than  it  should  be  by  a  small 
amount.  Convergence  studies  have  shown  that  as  the 
number  of  elements  on  the  body  surface  is  increased  the 
results  tend  toward  those  obtained  using  the  image  method 

The  inverse  of  the  normalwash  factor  matrix  [D]  need  not 
be  found  to  obtain  a  solution.  Direct  solution  by  triangu- 
larization  is  always  preferable  If  the  same  problem  is 
to  be  solved  repeatedly,  a  method  for  saving  the  triangu- 
larization  instructions  (these  are  numbers  which  may  be 
saved  in  the  lower  half  of  the  t riangulan-zed  [Dl  matrix) 
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Figure  1.  Variation  of  the  Total  Lift  Coefficient  with  Box  Number  for 
a  Wing  Tail  Combination  Oscillating  in  Plunge. 
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Figure  2.  Variation  of  the  Total  Lift  Coefficient  with  the  Vertical 
Distance  Between  Wing  and  Tail  for  Plunging  Motion. 


Figure  3*  Variation  of  the  Total  Lift  Coefficient  with  the  Vertical 
Distance  Between  Wing  and  Tail  for  Plunging  Motion, 


(b)  FUSELAGE  WITH  CAMBER,  WING  a=  0 


Figure  4.  Span  Load  Distribution  for  a  Wing-Body  Combination 
in  Steady  Flow;(a)  Body  Angle  of  Attack  =1,0  Radian, 
Wing  Angle  of  Attack  =  0.  0;  (b)  Body  with  Parabolic 
Camber,  Wing  at  Zero  Incidence. 


j>-a.-o-a-000i 

4  .Or  --O'  • 

J&"- — «-  1 

2.0.  "\N 


i  0 
-2.0 
-4.0 
-6-0 

)  6.0 1 

c. 


0.2  0.4  0.6 

y/* 


M  =  0 
*r=  1.72 


-o--  R'cAL  |  ' 

IMAGINARY 

* 

o  o  WING  AND  NACELLE|method 

•  «  WING  ALONE 


'i  4_0L'0'"-o--o-o-D-ooca 


0  02  0.4  06  06  10 

y/s 


Figure  7,  Distribution  of  Lift  Coefficient  for  a  Wing-Nacelle 
Combination-  (a)  Plunging;  (b)  Pitching. 
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Figure  8.  Comparison  of  Experimental  and  Calculated  Lifting 

Pressure  Coefficient  for  a  Wing-Nacelle  Combination 
in  Pvunge. 
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APPENDIX  A.  THE  KERNEL  FUNCTION  FOR  LIFTING  SURFACES 


The  velocity  normal  to  an  oscillating  surface,  W  =  U^Re^We1^),  is 

icofc 

related  to  the  pressure  difference,  P  =  q  Re(AC  e  w  ),  across  the  surface  by 

P 

the  integral  equation  (Ref.  16) 


w,  M)AC  (|,  r\,  C)d£  do¬ 
's 


(A.  1) 


where  x  or  £  is  the  streamwise  coordinate  and  <r  is  the  tangential  coordinate 
on  the  surface  (see  Figure  A-l),  w  is  the  frequency  of  oscillation,  M  is  the 
Mach  number,  and  subscript  r  denotes  the  downwash  or  receiving  point  and 
the  subscript  s  denotes  the  doublet  or  the  sending  point,  and 


xn  =  x  “  % 


(A.  2) 


Vn  =  V  -  1 


(A.  3) 


zQ  =  z  -  i 


(A.  4) 


The  symbol  J-  means  that  the  integral  in  Equation  (A.  1)  is  defined  in  the  sense 
of  its  "finite  part"  (Ref.  17). 


Rodemich  (Ref.  7  )  has  derived  an  expression  for  the  kernel  function 
for  a  nonplanar  surface  in  the  form 


K 


e-‘“x0/U.(K 


iTi/r' 


+  K 


* ,  4 
2T2/r  > 


(A.  5) 


where 


r 


<4 


2,1/2 

'0  ' 


T1  =  cos^r-  V 


(A.  6) 


(A.  7) 
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(A.  8) 


T2  =  (zQ  cos  Yr  -  yQ  sin  Vr)  (z0cos  Ys  -  yQ  sin  Ys) 


(N.B.,  the  definition  of  T'  in  Equations  (A.  5)  and  (A.  8)  differs  from  that  of  T~  m 

^  2 

Refs.  7  and  8  by  the  factor  r  .) 


(  8  ) 

and  Landahl'  '  has  simplified  the  forms  of  and  to  read 


where 


" i  ui  2  1/2 

Kj  =  Ij  +  [Mr  / R]  [•  e  /(I  +  Uj )  '  } 


(A.  9) 


K„  =  -3L 


ik1M2r2 

R2 


■iklul 


n  x  V/2 
(1  +  uL) 


Mr  j 

IT 


(1+u2) 


2  2 

P  r 

R2 


Mr  u 


+  2  + 


R 


(A.  10) 


e"iklUi 


u]L  =  (MR  -  Xq)/(3  r 


(A.  11) 


k,  =  wr  /U 
1 


(A.  12) 


P  =  (i-m2)1/2 


(A.  13) 


R 


,  2  Z  2.1/2 

=  (xQ  +  P  r  ) 


(A.  14) 


and 


-ikj  u 


Ij  (Uj,  kj)  = 


/  nx  2>3/ 2 

/Uj  (1+u  ) 


du 


(A,  15} 
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"2<V  fcl>  =  J 


(Hu2, 5/2 


Figure  A-l 


At  zero  frequency  the  planar  part  of  the  kernel  is 


K<S>  -  T  +  _  _J _ 

K‘  1  R  (H-u2)1'2 


where 


(l+u2)3/2 


=  1  -  Uj  (1  +  u2)~1/2 


so  that 


(A.  10) 


(A.  17) 


(A.  18a) 


(A.  18b) 


1  I  xn/R 


(A.  19) 


The  zero  frequency  value  of  the  nonplanar  part  of  the  kernel  is  easily 
found  from  the  planar  part  {Refs.  7  and  8), 


{A.  20) 


(A.  21) 


(A.  22a) 


=  -  2  -  (x0/R)  (2  +  P2r2/R2)  (A.  22b) 

Although  the  integrals  Ij  and  [Equations  (A.  15)  and  (A.  16)]  are  easily 
evaluated  at  zero  frequency,  no  explicit  solution  is  possible  for  nonzero 
frequency.  However,  the  symvnetry  properties  of  the  integrands  permit 
consideration  of  only  nonnegative  arguments  u^  since,  for  <  0,  we  note 
that 

Ij  (Uj,  kj)  =  2ReI1  (0,  k^)  -  R-eL^  (-u^  kj)  t  ilml^  (-Uj ,  kj.)  (A- 23) 

and 

I2  (u^,  kj)  =  2ReI2  {0,  kj )  -  Rel^f-u^,  kj)  +  iiml^f-Uj,  kj)  (A  24) 
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.1 


1 


'•  p . 


fr- 


i 


£ 


i 


The  integration  of  Ij  by  parts  once  gi 


ives 


-ik,  u  2-1/2  -i 

Ij  (u^  k^)  =  e  [1  -  u!  (1+ui)  ‘  lki  W  kl^ 


(A.  25) 


where 


r  ® 

ik,u.  /  2  - L /2  ,  “*klu 

0lul'  V  *  c  J  [‘-uU-u)  1« 


du 


(A.  26) 


Integrating  I2  by  parts  twice  leads  to 


■ik,  u. 


v.  lx.  (  7.1/2 

3I2(ult  kj)  =  e  {  (2+ikjUj )  [1  -  Ujd+Uj)  ] 


-  u^l+u2)-372  -  ikjl0 (Ul,  kt)  +  ^J0(u1, 


(A.  27) 


where 


W  ki’  =  c 


ik|U'  f~  „[1  -u(l+u2)-1/2]c"klU  du  (A.  28) 


'UJ 


2, -1/2 


The  integrals  IQ  and  JQ  can  be  evaluated  using  approximations  to  u(l+u  ) 
Laschka^1^^  has  obtained  an  extreme"  y  accurate  approximation  in  exponential 


form  for  u  £  0 


1  -  u(l+u2f1/Z  s  £  ane’ 

n=l 


ncu 


(A.  29) 


where  c  =  0.  372  and  the  an  are  given  in  Table  I.  The  maximum  error  of 
this  approximation  is  0.  135%.  The  integrals  JQ  and  JQ  for  u  ^  0  then 
become 


-ncu 


1 


W  ki>  *  E  Tfcj  <nc  •  ik>’ 


(A  30) 


,  n  c  +  k, 
n=  1  1 
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>A 


and 


11 

W  kr  2 

n=l 


?ncu 


a  e 
n 


1 


,2  2-  2.2 

(n  c  +  k^) 


2  2  ,2 

n  c  -  kj 


(A.  31) 


+ 


2  2 

ncu^(n  c 


+ 


2 

1 


ik[2nc  + 


,  2  2 
u^n  c 


+  #] 


1 

j 


The  desired  integrals  and  then  follow  from  Equations  (A.  25)  and  (A.  27) 
for  2  0,  and  from  Equations  (A.  23)-(A.  25)  and  (A.  27)  for  u^  <  0. 


Table  I  -  Coefficients  in  Laschka's  Approximation 
to  u(  1  +  u2)"1/2 


n 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


+  0.24186198 

2.7918027 
+  24. 991079 

111.59196 
+  271.43549 

305.  75288 
41.183630 
+  545.  98537 

644.78155 
+  328.72755 

64.279511 
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is  warranted.  The  value  of  the  arctangent  is  taken  in  the  range  (-it  »+  it), 
_2  _2  2 

When  y  +  z  >  e  ,  Equation  B.6b  is  well  behaved  and  approaches  zero 
as  z  approaches  zero.  Under  these  circumstances  it  is  convenient  to 
introduce  a  parameter  a  defined  by 


F  = 


2e 


2  L  _2  2 
y  +  z  -  e 


1  —  a 


_2 

z 


(B.  7) 


where 


4e 


ass 


£  — 


2ez 


2n-4 


,—Z  .  -2  2,2  .  ,  \_2  _2  2 

(y  +  z  -  e)  n  =  2  2n-l  \y  +  z  -  e  . 


(B.  8) 


3  2  2 

Eq.  (B.  8)  has  been  used  for  |  2  ez/ (y  +  z  -  e  )  j  £  0.  3  when 
_2  _2  2 

y  +  z  >  e  .  Then  the  planar  downwash  factor  becomes 


D 


(1)  _ 


rs 


Ax 
_ s_ 

8tt 


[(y2  -  z2)  Aj  +  yB,  4  Cj]F 


,2  .  _2 

,  (y-e)  +  z 

+  (~2  +  yAj  )  log  ~  Y  +  ZeAjUxJ.  9) 

(y+e)  +  z  ) 


The  incremental  nonplanar  oscillatory  downwash  factor  is  approximated  by 


D 


(2) 

rs 


Ax. 


8  TT 


/ 

-G 


_ P?(9) 

[(y-n)2  +  z2]2 


dn 


(B.  10) 
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where  F* ^ ( T1 )  is  another  parabolic  approximation 


Pj  (n) 


B2n  +  c2 


(B.  1 1  a) 


{K2  exp  [ -iw(x-n  tan  X  )/U  ]  -  K 


(s) 


}t. 


(B. lib) 


Letting  P2(-e),  P2(0),  and  P2(e)  denote  the  inboard,  center,  and  outboard 
values  of  P2(*l),  respectively,  we  have 


A2  =  [P2(-e)  -  2P2(0)  +  P2(e)]/2e^ 


(B. 12) 


B2  =  [P2(e)  -  P2(-e)]/2e 


(B. 13) 


c2  =  P2(0) 


(B. 14) 


The  nonplanar  downwash  factor  is  then  evaluated  to  be 
Ax 


D 


2rs 


1 


(y2j-z^)A2  i  yB2  +  C2 


(B. 15) 


(y+e) 


~ — f[(y2+z^)y  +  (y2-z2')e]A2  +  (y2+z2+ye)B2  +  (y+eJC^ 
+  z  '  ' 


(y-e) 


jT - 2  (  C  (y^+z21)  y  -  (y2--z2)e]A2  +  (y2+z2-ye)B2  +  (y-e)  C2 

+  z  ' 


Eq.  (B.  15)  tends  to  lose  significance  for  small  values  of  z.  Introducing 
a  into  Eq.  (B.  15)  thru  Eq.  (B.  7)  leads  to 


Ax„ 


_  ,-2.-2  2. 

8ir(y  +  z  -e  ) 


(  2(y“+z2+e2)  (e^A2+C2)  +  4ye2  B2 
(  [  (y-e)2  +  z2]  [  (y-e)2  +  z2  ] 


-  “7  [  (y2+z2)  A?  +  yB 


2  +  C2  ^ 


(B. 16) 


The  simplification  of  Eq.  (B.  15)  via  Eq.  (B.  7)  is  somewhat  tedious  but 
results  in  the  more  accurate  form  above  in  which  a  is  again  given  by 
Eq.  (B.  8).  Eq.  (B.  16)  has  been  used  in  general  except  when 
|  (y- +  z  -e  )  /  2ez  |  S  0  .  1  in  which  case  Eq.  (B.  1 5)  is  used. 

2  2  2 

When  y  +  z  <  e  then  the  expression  for  F  is  not  well  behaved  but 
approaches  infinity  as  z— - 0.  The'specific  behavior  is 

lim  „  .1-1  ,  _2  _2  2 

F  —  ir/|z|  ,  for  y  +  z  <  e 


Under  these  circumstances  it  can  be  shown  that  the  term 


lim 
z  •— 0 


F  jcj  4  yB,  -1  (y2-z2)A1  +—  (C2  +  yB2  +  (y2+z?')A2)| 


—  0 


even  as  F  — ■  a>  ,  Numerically,  of  course,  this  expression  will  lose  accuracy 

X“  ^ 

at  some  value  of  z.  An  example  calculation  has  shown  that  for  - =  0,  5, 

k  =  6.  0  and  e  =  1/1 20, the  gap,  z/e,  at  which  numerical  difficulties  arose  was 
0.  0025.  Since  e  is  small  z  =  e  (.  0025)  is  very  small  and  can  be  assumed  to 
be  zero;  i.  e.  ,  the  planar  case. 
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APPENDIX  C.  THE  STEADY  NONPLANAR  DOWNW ASH  FACTOR 


The  steady  nonplanar  downwash  factor  is 

,  .  e  1  K*S*T  \ 

D<ss>  =  <a%/8„)  r  ^--V-  +  -jy  «  <cl> 

The  downwash  factor  has  been  given  by  Hedman  (Ref.  2,  App.  C)  and  is 
equivalent  to  Equation  (C.  1),  although  it  is  derived  from  fundamental  vortex 
considerations  rather  than  evaluating  the  integral  in  Equation  (C.  1).  We 
find  it  convenient  to  rederive  .Hedman' s  results  in  order  to  facilitate  the 
programming  of  the  equations. 

The  vector  geometry  for  a  finite -length,  constant-strength  vortex  is 
shown  in  Figure  C-l.  The  vector  form  of  the  Biot-Savart  Law  for  this 
segment  is 


V 


r  x  d 

4ir  dZ 


(cos  0 


cos  <P  ) 


(C  2) 


where 


d 


(  r/  r)  r^  cos  0 


(C  3) 


cos  0  =  (SjVRj.)--  •  (r/r) 


(C.4) 
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and 


cos  cp  =  (R2/R2)  .  (T/r)  (C.5) 

and  where  an  arrow  indicates  a  vector.  The  normalwash  at  the  receiving 
point,  where  the  dihedral  angle  is  y  ,  is 

w  =  V  .  (j  sin  Vr  -  k  cos  y  )  (C.  6) 


where 


v  =  jvy  +  kVz  (C.7) 

The  general  result  desired  for  a  horseshoe  vortex  system  applies  the 
foregoing  equations  to  the  bound  and  two  trailing  vortex  components  of 
the  horseshoe  system. 

The  physical  geometry  of  a  horseshoe  vortex  system  is  shown  in 
Figure  C.  2(a);  the  geometry  modified  by  the  Prandt-Glauert  Transformation 
is  shown  in  Figure  C.  2(b).  The  transformations  are 

x0  =  xQ/p  (C.  8) 


and 


tan  A  =  (tan  y  )  /P 


(C.  9) 


The  Biot-Savart  Law  is  to  be  applied  using  the  geometry  of  Figure  C.  2(b). 
Let  the  origin  of  coordinates  be  the  center  of  the  bound  vortex.  Then  the 
distance  to  the  receiving  point  (Xq,  y^,  Zq)  is 


R  =  i  XQ  +  jyQ  +  kzQ 


(C.  10) 
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and  the  distances  from  the  inboard  and  outboard  ends  of  the  bound  vortex 
are 


r . 
1 


-  e  (i  tan  A  +  j  cos  y  +  k  sin  y  ) 

s  s 


(C.  11) 


and 


e  (i  tan  A  +  j  cos  y  +  k  sin  V  ) 

s  s 


(C.  12) 


respectively.  The  distance  to  the  receiving  point  from  the  inboard  end 
is 


R.  =  R  -  r. 
i  i 


(C. 13a) 


=  i  (XQ  +  e  tan  A  )  +  j  (yQ  +  e  cos  y&)  +  k(zQ  +  e  sinyg) 

(C. 13b) 

*  T  Rbc  +  TRiy  +  “Sj  «=•  13=) 


from  which 


R.  = 


(R2  +  R2  +  R2  )1/2 

ix  ly  iz  ' 


(C.  14) 


The  distance  to  the  receiving  point  from  the  outboard  end  is 


R  =  R  -  r 
o  o 


(C.  15a) 


i(XQ  -  e  tan  A)  +  j(yQ  -  e  cos  Yg)  +  k(zQ  -  e  sin  y  )  (C.  15b) 


iR  +  jR  +  kR 
ox  J  oy  oz 


(C. 15c) 
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from  which 


R  =  (R2  +  R2  +  R2  )1/2  (C.  16) 

o  ox  oy  oz  ' 

The  vorticity  vectors  for  the  three  components  of  the  horseshoe  vortex 
complete  the  basic  data  required  for  calculation  of  the  downwash  factors, 

r,  -  T(i  sin  A+  j  cosAcos  V  +  k  cos  A  sin  7  )  (C.  17) 

b  s  s 

r.  =  r(-T)  (c.  18) 

ro  =  r(+T)  (c.  19) 


We  begin  the  final  calculations  by  finding  the  cosines  of  0  and  <P  and 
the  normal  distance  vectors  d  for  the  bound,  inboard,  and  outboard  vortices. 


For  the  bound  vortex, 

cos  eb  =  (R./Ri)  .  (rb/r)  (C.20a) 

=  (R.  sinA  +  R.  cos  A  cos  y  +  R.  cosAsin  y  )  (C.  20b) 

R.  '  ix  xy  s  iz  s 

i  1 

cos  <Ph  =  (R0/RQ)  .  (rb/T)  (C.  21a) 

=  jr  (R  sin  A  +  R  cosAcos  Y  +  R  cosA  sin  V  )  (C.  21b) 

R  '  ox  oy  s  oz  s 

o 

d“b  =  R.  -  (rb/r)R.  cos  0b  (C.  22a) 


=  i(R.  -  R.  cos  0,  sin  A) 
'  ix  l  b  ' 
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;  i 

i  +  j(R;,,  “  R.  cos  0,  cos  A  cos  y  ) 

iy  i  b  s 


l 


l 


I  | 
t  f 


t:  ; 


¥• 


=  idbx  +  Jdby  +  kdbz 


For  the  inboard  vortex. 


r. os  0.  =  +1 


cos  <?. 
1 


(Rj/Rj)  .  (r./r) 


-  -  R.  /R. 

IX  x 


\  =  jRiy  +  kR 


1Z 


Finally,  for  the  outboard  vortex, 


cos  0 


<VRo>  •  l?o/r) 


=  R  /R 

OX  0 


cos  <P  =  -1 

o 


d  =  jR  +  kR 
o  J  oy  oz 


We  next  evaluate  the  velocity  induced  by  each  of  the  three  vortices.  The 
bound  vortex  induces  the  velocity 


Fb  X  dh 
4*d  2 


(cos  0b  -  cos  <P  ) 
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(C.  22b) 

(C. 22c) 

(C.  23) 

(C.  24a) 

(C.  24b) 

(C.  25) 

(C.  26a) 
(C. 26b) 
(C.  27) 

(C.  28) 

(C. 29a) 


=  — — 2~  (cos  Gb  -  cos  <»b)  [i  (dbz  cos  Acos  -  d.  cos  Asin  yg) 

4ndb 

+  j(^jjX  cos  Asin  7s  -  dbz  sin  A  )  (C.  29b) 

-r  "k(db  sin  A  -  dbx  cos  Acos  yg)  ] 

=  r(Tvbx  +  TVj  ,  +  ■£vbz)  (C.  29c) 


(C.  30) 


(C.  31a) 


(C. 31b) 

(C.  31c) 


(C. 32a) 


(C.  32b) 


(C. 32c) 
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The  total  velocity  induced  by  the  horseshoe  vortex  is 


V  =  V,  +  V.  +  V  (C.  33a) 

bio 

=  rrT(v.  +  v.  +  v  )  +  T(v,  +  v.  +  v  ) 

L  bx  ix  ox  •'  by  ly  oy 


+ 

k(V,  +  V.  +  V  )  ] 
bz  iz  oz  ‘ 

(C  33b) 

- 

r(Tv  +  Jv  +  "kv  ) 

'  x  J  y  z 

(C. 33c) 

The  normalwash  velocity  is  then  given  by 

w  = 

V  .  (j  sin  -  k  cos  y^) 

(C.  34a) 

= 

r(V  sin  y  -  V  cos  y  ) 

y  X  Z»  Jl 

(C.  34b) 

which  leads,  finally, 

to  the  steady  nonplanar  downwash  factor 

D(S) 

=  Ax  w  /2r 

(C.  35a) 

rs 

s  r 

=  (V  sin  y  -  V  cos  y  )  Ax  / 2 
'  y  r  z  r  s 

(C.  35b) 

A  numerical  difficulty  arises  when  the  receiving  point  is  on  an  extension 

of  the  bound  vortex  axis.  The  limit  of  as  d^  approaches  zero  is  also 

zero.  Numerical  accuracy  can  be  preserved  by  utilizing  a  series  expansion 

of  (cos  0,  -  cos  <P,  ).  For  small  0,  and  <P, 
bb  b  b 


cos  0,  -  cos  <P. 

b  b 


°b>'2 


(C.  36) 
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where 


In  either  case 


(cos  0. 
b 


cos  9^)  /  d 


2 

b 


(C.  41) 


The  program  uses  the  apj.  “  iximation  if  max  (  |cos  0^  |.  |  cos  9^  |)  >  0,999- 

If  0^  «  0  and  9^  ~  w,  the  receiving  point  is  close  to  the  bound  vortex  This 

case  is  only  of  interest  in  induced  drag  calculations.  Then  the  bound  vortex 

2 

induces  no  velocity  on  itself.  The  program  sets  (cos  0^  -  cos  9^)  /  d^  to 
zero  if  cos  0^  cos  0  and  max  (  jcos  0^  |  ,  |  cos  9^  | )  >  0-  999, 


APPENDIX  D.  THE  SUBSTANTIAL  DERIVATIVE  AND  INTEGRATION  MATRICES 


D.  1  The  First  Alternative 

Alternative  No.  1  for  the  substantial  derivative  matrix  [W]  and  the 
integration  matrix  [B]  assumes  that  the  deflections  of  a  given  spanwise 
division  (aerodynamic  strip)  can  be  determined  by  the  following  degrees  of 
freedom. 

h  -  plunging  of  strip 

a  -  pitching  of  strip 

p  -  rotation  of  primary  control  surface 

6  -  rotation  of  secondary  control  surface 

P  -  plunging  of  primary  control  surface 

6  -  plunging  of  secondary  control  surface 

To  describe  the  format  of  the  matrices  involved,  we  assume  that  a  spanwise 
division  (strip)  is  idealized  into  10  chordwise  boxes  as  shown  in  Figure  D-l. 
Any  number  of  boxes  may  be  chosen;  10  are  used  only  as  an  example.  The 
local  deflection  mode  shapes  corresponding  to  the  local  degrees  of  freedom 
h,  a,  p,  6,  p,  and  6  are  also  shown  in  Figure  D-l.  The  [W]  and  [B] 
matrices  are  constructed  from  the  following  equations. 


ti.  =  h  +  (x.  -  x. )  a 
J  J  1 


X.  <  x„ 

J  2 


(D.l-l) 


(D.l-2) 


h+  (x  -  x  )«+  (x  -  x  )p  +  pj 

J  J  (  x  2  <  xj  <  x4 


h  +  (x.  -  x  )v  +  (x.  -  x,)P  +  p 
J  i  J 

+  (Xj  -  x4)6  +  6 


=  a  +  p  +  6 


,Xj  >x4 


(D-  1-3) 


(D.  1  -4) 
(D.l-5) 


(D.l-6) 


?-  PLUNGING 


PITCHING 


III 


i.  -  Typical  distance  to  1/4  chord  of  "box"  f  ° 

j 

Xj  -  Typical  distance  to  3/4  chord  of  "box" 

-  Distance  to  reference  axis 

x^  -  Distance  to  leading  edge  of  primary  control  surface 
x^  -  Distance  to  hinge  line  of  primary  control  surface 

-  Distance  to  leading  edge  of  secondary  control  surface 
Xg  -  Distance  to  hinge  line  of  secondary  control  surface 

Figure  D-T.  Strip  Geometry  and  Local  Deflection  Shapes  for 
Alternative  1 
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The  format  of  the  [W]  matrix  which  results  from  the  above  equations 
for  the  example  of  Figure  D-l  is  shown  in  Table  D-l.  The  format  of  the 
[B]  matrix  is  not  shown,  but  can  be  determined  by  substituting  £.  for  x. 
in  the  [W]  matrix  and  writing 


[B] 


(D.  1 


where  [Aj  is  the  diagonal  matrix  of  box  areas.  The  expression  q  [B],  then 
represents  lifts  due  to  h,  (3,  6,  and  moments  (about  x^,  x^,  x,.)  due  to  a,  p, 
6. 

The  present  version  of  H7WC  generates  AIC's  only  for  lifting  surfaces; 
fuselage  elements  or  bays  are  ignored.  (See  Part  I,  Vol  II)  . 
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TABLE  D-l.  MATRIX  [W]  FOR  TYPICAL  STRIP 
ALTERNATIVE  NO.  1 


The  Second  Alternative 


Alternative  No.  2  which  includes  parabolic  camber  is  shown  in 
Figure  D-2.  The  Lagrangian  interpolation  formulas  are  used  to  describe 
the  deflections  of  the  strip.  The  applicable  formulas  are: 

(x.  -  x,)(x;  -  x,)  (x.  -  x.  )(x.  -  x,) 

v,  =  J  -  ,i  $  u  ±  1  1  3 _ _ _  b 

j  (xL  -  x2)(x1  -  x3)  1  (x2  -  x1)(x2  -  x3)  2 


(x.  -  x.  )(x.  -  x_)  )  x.  <  x_  when  control  surface 

+  t-^ — - — — — rh,,|  ^  -is  present; 

'x3  X1  x3  x2  all  x.  with  r.o  control  surface. 

)  J 


(D.2-1) 


dh.  /(x.  -  x2)+(x.  -  x3) 


_ L  -  _J _ £_ 

dx  \(x  -  x2)(x 


c.  -  x,)  \  / (x.  -  x.  )+(x.  -  x 

J - 2-  h,  +  Ir-J - -  1 - 

1  -  x3)  /  1  \(x2  -  xx)(x2  -  x 


"  X3} )  h 
■  *3)  /  2 


(X.  -  w1)+(x.  -  x2) 
(x3  -  x1)(x3  -  x2) 


))  X.  <  X,  V 

h3’  |  1 

j  all  Xj  wi 


x.  <  x3  when  control  surface 
is  present  ; 

all  x.  with  no  control  surface. 


(D.  2-2) 


When  the  control  surfaces  are  present,  we  assume  that  the  deflections  aft 
of  x3  can  be  determined  by  linear  extrapolation.  We  then  have 


j(y.3  -  x2)(x.  -  x3)  \  ^  ( (x3  '  '  x3)  \  i 

hj  "  y(x3  -  x][)(x1  -  x2)j  hl  +  ^(x2  -  x1)(x2  -  x3)  J  h2 

/  /{x  -  x  H(x  -  x?)  \  \ 

+  I1*  <*j  -  *3>  )  hy  >  *3  (D.  2-3) 


„  h  |X3  -  *2»  .  .  <x3  -XI> 

1  (x1  -  x2)(Xl  -  x3)  +  1  2  (x2  -  Xl)(x2  -  x3) 


(x3  -  x1)+(x3  -  x2) 

3  (x3  -  x^^  -  x2)  '  xj  >  x3 


(D.  2-4) 
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Figure  D-2.  Strip  Geometry  and  Local  Deflection  Shapes  for 
Alternative  2. 


The  incremental  deflections  due  to  the  control  surface  degrees  of  freedom 
are  given  by  the  equations  below. 


Vj 

=  (*j  -  x5)  p, 

x.  >  x . 

J  4 

(D.  2-5) 

3(A,h.) 

3x 

=  p. 

x.  >  x . 

J  4 

(D.  2-6) 

Ah. 

2  J 

=  (Xj  -  X?)  6, 

X.  >  X, 

J  7 

(D.  2-7) 

3(A2h.) 

3x 

=  6, 

x.  >x7 

(D.  2-8) 

A3hj 

=  P> 

Xj>X4 

(D.  2-9) 

A4"j 

=  6, 

Xj>X7 

(D.2-10) 

The  [W]  and  [B]  matrices  can  now  be  constructed  from  Equations(D.  2-1) 
through  (D.2-10).  Figure  D-2  indicates  in  a  qualitative  manner  the  local 
deflection  mode  shapes  corresponding  to  the  degrees  of  freedom  shown. 
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